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Abstract 


Simulations of nonlinear conservation laws that admit discontinuous solutions are 
typically restricted to discretizations of equations that are explicitly written in diver- 
gence form. This restriction is, however, unnecessary. Herein, linear combinations 
of divergence and product rule forms that have been discretized using diagonal- 
norm skew-symmetric summation-by-parts (SBP) operators, are shown to satisfy 
the sufficient conditions of the Lax-Wendroff theorem and thus are appropriate for 
simulations of discontinuous physical phenomena. Furthermore, special treatments 
are not required at the points that are near physical boundaries (i.e. , discrete con- 
servation is achieved throughout the entire computational domain, including the 
boundaries). Examples are presented of a fourth-order, SBP finite-difference op- 
erator with second-order boundary closures. Sixth- and eighth-order constructions 
are derived, and included in E. Narrow-stencil difference operators for linear vis- 
cous terms are also derived; these guarantee the conservative form of the combined 
operator. 
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1 Introduction 

The nonlinear product terms in the continuous Euler equations can be expressed 
in many forms: 1) conservative, ^r{uiUj)\ 2) primitive, 3) skew-symmetric, 

l£]( u i u j) + 4) rotationa1 ’ u i(t| “ fe) “ £](l u i u iY’ as wel1 as o^- 

ers. Although they are equivalent for smooth flows at the continuous level, each 
form can exhibit profoundly different discrete accuracy, conservation, and nonlin- 
ear stability properties. Growing evidence suggests that the nonconservative forms, 
particularly the skew-symmetric form, preserve important invariants of the contin- 
uous equations (e.g., kinetic energy or entropy) and deliver enhanced accuracy and 
robustness. This makes the use of split forms an attractive alternative to conven- 
tional approaches based on the conservation form of the equations. Evidence that 
supports this assertion is available in the following references [1—16]. 

Despite many favorable properties, alternative, nonconservative forms of the 
Euler equations have historically been viewed as inappropriate for simulations that 
involve discontinuities, such as supersonic and hypersonic flows. The basis of this 
misconception likely stems from the Lax-Wendroff theorem [17] which stipulates 
that the divergence form of the governing equations must be discretized with a 
conservative operator, a condition that alternative discrete forms apparently lack. 1 

While alternative forms are generally not explicitly conservative, many can be 
manipulated into an equivalent, consistent, and conservative form, thus allowing 
the Lax-Wendroff theorem to be used to guarantee that the convergent captured 
discontinuities are weak solutions of the governing equations. For example, Ducros 
et al. [9] demonstrated that the skew-symmetric form of the convective terms in the 
Euler equations are discretely conservative (at least for periodic fourth- and sixth- 
order centered operators), while Jameson [18] used the skew-symmetric form of the 
Burgers equation and developed discretely conservative second order operators that 

1 Note that Gerritsen and Olssen [6] and, later, Yee and Vinokur [10] showed that correct shock 
simulations could be achieved using the skew-symmetric form, even before a proof of discrete 
conservation was identified. 
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are capable of capturing shocks. Recently Pirozzoli [15, 16] showed that all periodic 
centered finite-difference operators can be expressed in a conservative form (even on 
curvilinear grids) if an appropriate splitting of the equations is used. 

Several open issues must be resolved before high-fidelity simulations that use 
alternative forms of the equations become routine. One critical requirement is the 
need for high-order boundary closures and boundary conditions that maintain the 
desirable properties of the interior operators. Another is the need for a systematic 
methodology to extend alternative single-domain operators (e.g., a skew-symmetric 
discretization) to multiple domains, thereby extending the generality of any alter- 
native approach. 

A principle contribution of this work is the recognition that the boundary clo- 
sures for alternative operators follow immediately if finite-domain summation-by- 
parts (SBP) discrete operators are used to derive near-wall boundary closures. For 
example, boundary closures that preserve the skew-symmetric form of the Euler 
equations up to and including the boundary conditions are achieved by appropriately 
splitting the equations and using SBP operators for all derivatives. We demonstrate 
this point by developing three sets of fully conservative, alternative finite-difference 
operators in this study: (2-4-2), (3-6-3) and (4-8-4). 2 

Another contribution of this work follows from addressing the broader question, 
“Which discretizations of general conservation laws are discretely conservative?” 
To this end, we conjecture and support with many examples that any combina- 
tion of the conservative and primitive conservation law operators, otD(yw ) + (1 — 
a) [vT>w + wT>v), when discretized using any , diagonal-norm, skew-symmetric, SBP 
discretization operator, can be expressed in a conservative form after the appropriate 
discrete manipulations have been made. The generality of this conjecture allows for 
quite arbitrary splittings of the equations while maintaining discrete conservation. 

A final contribution of this work is the recognition that all single-domain al- 
ternative operators satisfy all of the conditions that are necessary for extension 
to general, multiblock discretizations, provided that interface coupling terms that 
maintain stability, accuracy, and conservation are used [21]. 

Another critical research element is the need for dissipation operators that are 
complimentary with the skew-symmetric form, but maintain the design-order accu- 
racy of the finite-domain operator. Rather than include this topic herein, we address 
it in a companion paper. Because all of the discrete operators that are derived in this 
paper are nondissipative and have no inherent mechanism to suppress oscillations, 
we refrain from testing problems that admit strong shocks. As such, we validate the 
accuracy of the new class of operators by using the viscous Burgers equation. 

The layout of the paper is as follows. Section 2 defines SBP operators, briefly 
summarizes the Lax-Wendroff theorem, and presents a derivation of the discretely 
conservative split-form finite- difference approximations for conservation laws. Ex- 
amples of a complete fourth-order operator are presented in section 3. Section 3 also 

2 The ( p-2p-p ) nomenclature denotes that boundary-interior-boundary stencils are p-, 2 p-, p- 
order accurate, respectively. The resulting operators are (p+l)-order accurate for hyperbolic equa- 
tions (Euler) and incompletely parabolic equations (Navier-Stokes) and (p+2)-order accurate for 
parabolic equations (viscous Burgers) [19,20]. 
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contains a discussion on the relationship between the SBP operators and split-form 
conservation, as well as a discussion of how alternative operators can be extended 
to complex geometries. In section 4, several split-form operators are applied to the 
viscous Burgers equation. Numerical validation of accuracy of an alternative form 
is given in section 5, with the use of both single and multiple domains. For brevity, 
the (3-6-3) and (4-8-4) operators appended verbatim in E. 

2 Definitions 

2.1 Summation-By-Parts Operators 

Discretize the physical domain by using N equidistant points that are distributed 
over the interval [0, 1] as 

x = (xi,x 2 , ■ ■ • ,x N -i,x N ) , x i = (i - 1) , i = l,2,...,JV 

and define a general class of SBP derivative operators T> sbp that satisfy the following 
matrix properties: 


V sbp = V- 1 Q, Q + Q t = B = diag(-l, (),••■ ,0,1) 

P = diag (p( U )r •• ,P(jv,at)), C T PC > 0, C/° 

where P(i.i) • ■ ■ P(n,n) are the diagonal elements of the matrix and £ is an arbitrary 
vector. The SBP operators defined in eq. 2.1 utilize a diagonal norm V and are, 
therefore, a subset of the more general SBP form (see ref. [21] or ref. [22]). All are 
telescoping operators and are conservative in the V norm. (See A for a proof that 
all operators that can be expressed in the form of eq. 2.1 are conservative.) 


2.2 Complementary Grids 

The discrete derivative operator approximating h x ( h sufficiently smooth) is repre- 
sented herein as a combination of matrix operations that transfer data between two 
complementary meshes. This nomenclature automatically leads to discrete conser- 
vation. The two complementary sets of grid points (or meshes) are defined on the 
finite interval 0 < x < 1, differ in dimension by one and are expressed by using the 
vectors 

x = [aq,X2,--- , xn-1, xn] T , x = [xq,x\,--- , xn-i, xn] T 


with x\ = xq = 0 and xjv = xn = 1. Herein, the points x are referred to as 
the “solution points,” while the points x are referred to as the “flux points.” (The 
overbar nomenclature is reserved herein for those quantities that are defined at the 
flux points x.) 

Define the finite-domain conservative difference operator in matrix notation 


as 


6h 

5x 


P _1 Ah 
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where 


/ —1 1 0 
0-11 


0 0 0 

V 0 0 0 


0 0 0 \ 

0 0 0 

0 0 

-1 1 0 

0 -11/ 


(2.2) 


is a rectangular (7V+1 x N) matrix. The flux vector h constructed using appropriate 
interpolation operators on the flux points x, is differenced back onto the solution 
points by A. The matrix V accounts for the appropriate grid spacings through- 
out the domain. Note that discrete conservation follows immediately because the 
columns of A sum to [— 1, 0, • • • ,0,1]. Further details on the definitions of the dual 
grids, as well as implementation details can be found elsewhere [24]. 


2.3 The Lax-Wendroff Theorem 


Consider a system of conservation laws of the form 


m + nzr = 0’ * > °> o < x < l, 

v ( x , 0 ) = vq(x) , v(0,t) = bc(t) 


(2.3) 


where v = v(x, t ) is the continuous solution vector, / = f(v) is a nonlinear flux vec- 
tor, and vq(x) is a bounded piecewise continuous function in L 2 . Although smooth 
solutions to eq. (2.3) do exist, no smoothness can be guaranteed due to the non- 
linearity of /. Thus, we seek a weak solution to eq. (2.3) such that the integral 
relation 


T 1 

■ 1 

T 

- rp 

/ / [i>tv + ip x f(v)] dxdt. = 

/ i fvdx 

+ 

[ iff(v)dt 

J d 
0 0 

J 

.0 

0 

J 

.0 


is satisfied for all smooth test functions i p(x,t). 

Choose a vector-valued function g(u) that is both a function of 21 compactly 
supported vector arguments and related to the flux / with the sole requirement 
that g ( v , ■ ■ ■ ,v) = f(v). Furthermore, define the discrete fluxes g to be of the form 


g(x + Sx/2) = g(u-i+i, ■ ■ ■ , u t ), g(x - Sx/2) = g(u-i, • • • , u t - 1 ) 


with suitably modified one-sided expressions near the physical boundaries. 

With these definitions, we state here without proof the Lax-Wendroff theorem 
[17] (originally derived for an infinite domain): 

Theorem 2.1. Consider the conservative discrete analogue of eq. (2.3): 

u(x,t + 8t)-u(x,t) g{x+ ^f) - g{x - f) _ (c? ^ 

St + Sx 

If the discrete solution u(x,t) converges boundedly almost everywhere to u(x,t), then 
u(x, t) is a weak solution to the continuous equation. 

Finite-domain proofs of the Lax-Wendroff theorem are given in the context of 
SBP operators in refs. [21,23]. 
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2.4 Split-Operator Consistency Requirements for Lax-Wendroff 

Consider the flux / that is defined in eq. (2.3), which satisfies the property f(u) = 
v(u)w(u). The general split flux which is composed of a linear combination of the 
divergence and the product-rule components, can be written as 

ut + af(u) x + (1 - a) [v(u)w(u) x + w(u)v(u) x ] = 0 (2.6) 

Discretizing eq. (2.6) with an SBP spatial operator results in the following semi dis- 
crete equation: 

u t + aV-'QWv + (1 - a) [VV^Qw + WV^Qv) = 0 

V = [v(ui),v(u 2 ), ■ ■ ■ ,v(u N )] T , V = diag(v) (2.7) 

w = [w(ui),w(u 2 ), ■ ■ .,w(u n )] t , W = diag(w) 

Equation 2.7 does not explicitly lead to a discretely conservative flux form i 2 . 
Furthermore, if the resulting discrete flux can be manipulated into a conservation 
flux that satisfies g(u, ■ ■ ■ , u) = f(u), then g is not immediately consistent with the 
physical flux / = vw that satisfies the Rankine-Hugoniot relations. Thus, if the 
Lax-Wendroff theorem is to be used to guarantee convergence to the weak solution 
of the continuous equations, then the following conditions must be met: 

1. The discrete spatial operator aV l QWw + (1 — a) (ViP^Qw + WV~ l Qv) 
must be equivalent to a telescoping form || . 

2. The discrete flux g is Lax-consistent with the physical flux /, where / is the 
flux that satisfies the Rankine-Hugoniot relations S[u] — [ f } = 0. (The bracket 
notation [] denotes the jump in the function.) 

3 A Discretely Conservative, Finite-Domain, Split Form 

3.1 A Fourth-Order Example 

We begin by demonstrating that the fourth-order finite-domain split-form operator 
is conservative for any value of the parameter a if a diagonal-norm, finite-domain 
SBP operator is used for discrete differentiation. 

The conservation part 'P~ 1 <2(!Fv) at an interior point i (s < i < N + 1 — s) 
is obviously conservative, and therefore does not need to be further manipulated. 
(The number of closure terms that are needed at each boundary is given by the 
constant s.) The product-rule terms VV~ l Qw + WV v Qy, which are discretized 
by using the conventional centered finite-difference at an interior point i. are 

(VV~ l Qw + WV~ l Qv)i = ^ ( Wi -2 - 8wj_i + 8w i+ i - w i+2 ) 

w t , , (3-1) 

+ — (Vi- 2 - 8Vi-! + 8v i+1 - V i+2 ) 
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and can be expressed in the following conservative form: 


( Vi ~ 2 \ 

T 

/o 0 1 0 0 \ 


/ VJi- 2 \ 

Vi - 1 


00-800 


Wi - 1 

Vi 


1-8 0 8-1 


Wi 

Vi+1 


0 0 8 0 0 


W i+ l 

\Vl+2/ 


\0 0 -1 0 0 / 


\Wi+ 2) 


fi = ~ 
ji 12 


f % 1 12 


eq. 3.2 are given as 





fVi- 2 \ 

T 

(0 

0 

0 

0 

0 ^ 


/ Wi- 2 \ 


Vi - 1 


0 

0 

0 

-1 

0 


Wi - 1 


Vi 


0 

0 

0 

8 

-1 


Wi 


Vi+1 


0 

-1 

8 

0 

0 


Wi+l 


\Vi+2j 


\0 

0 

-1 

0 

0) 


\Wi+ 2) 


/ Vi- 2 \ 

T 

(0 

0 

-1 

0 

0\ 


( Wi- 2 


Vi-1 


0 

0 

8 

-1 

0 


Wi- 1 

Vi 


-1 

8 

0 

0 

0 


Wi 

Vi+1 


0 

-1 

0 

0 

0 


Wi+l 

\Vi+ 2) 


\0 

0 

0 

0 

0) 


\Wi+ 2 

) 


(3.2) 


(3.3) 


Clearly, the coefficients in the flux matrix simply shift down and to the right for 
each subsequent flux fi (s < i < N + 1 — s). The boundary fluxes are more 
elaborate but telescope in the same manner. The boundary fluxes are completely 
determined by the form of Q near the boundaries for both the conservation and 
product-rule forms. 

The discrete split operator that is defined in eq. 2.7 can be written in con- 
servative form by using a linear combination of the conservative divergence and 
product-rule fluxes as 


u t + V x Af = 0, f = af c + (1 — a) f e 
The divergence form flux is given by 


/ 


fr = 


W\V\ 


gg (48uqui + 59^2^2 - 8 W 3 V 3 - 3W4V4) 

^ (-lluqwi + 59w 2 V 2 + 51w 3 v 3 - 3 W4V4) 
gg (-3uqui + 5 IW 3 U 3 -I- 56rc4U 4 - 8w 5 v 5 ) 

J2 (~w 3 v 3 + 7WiV4 + 7 w 5 v 5 - w e v 6 ) 

+2 (-Wi- lVi_l + 7 WiVi + 7w i+ lV i+ l - W i+ 2V i+ 2) 

J2 (~ w N-5VN-5 + 7WN-AVN-4 + 7wn~3VN~3 ~ WN-2 VN-2 ) 
(— 8un-4Vn-4 + 56u N - 3 v N - 3 + 51un-2Vn-2 - 3unvn) 
± (~3w N - 3 vjv-3 + 51WJV-2VN-2 + 59w N -iV N -l - Uw N v N ) 
^ (Swjv-^jv -1 - 8 w N - 2 v N -2 + 59u>Ar_iVAT — 1 + 48 w N v N ) 

W N V N 


(3.4) 


(3.5) 
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while the product-rule flux is 


W\V\ 


1 


— ( 59 u 2 Ui — 8113V1 — 3-U4U1 + 59 iiiU 2 — 8U1V3 — 3U1V4) 
96 

— (-8U3V1 - 3 « 4 Ui + 59 U3V2 - 8U1V3 + 59 tt 2 U 3 - 3 «ii> 4 ) 
9 p 

— (— 3114V1 + 59U4U3 — 8U5V3 — 3 rtiU 4 + 59-U3U4 — 8U3V5) 
96 

— ( VJ5V3 - W 3 V 5 + 8w 5 V4 + 8W4V5 - W 6 V4 - W6V4) 


l = 


12 


— ( - Wi-iv i+ i - w i+ iVi-i + 8 WiV i+1 
+ 8Wi + Wi - WiVi + 2 - Wi+2Vi ) 


( — WN-5VN-5 - WN-3VN-5 + 8WN-4VN-3 
+ 8WN-3VN-4 — WN-3VN-I - WN-1VN-3 ) 

1 . 

— ( — 8 UN- 2 VN -4 + 59iiAT_2UAr_3 — 3 unV ^-3 
96 

— 8UN-4VN-2 + 59UN-3VN-2 - 3UN-3VN ) 

1 . 

— ( — 3UNVN-3 + 59un-1VN-2 — 8UNVN-2 
96 

+ 59UN-2VN-1 - 3ujy-3VN — 8 UN- 2 VN ) 

1 . 

— ( - 3unvn-3 - 8unvn-2 + 59unvn-i 
96 

— 3un-3Vn — 8un-2Vn + 59un-ivn ) 
W N VN 


(3.6) 


We also have constructed operators for 2p = 6 and 2p = 8 with boundary closure 
blocks of dimension s = 7 and s = 9, respectively. The interior flux forms for the 
sixth- and eight-order schemes are presented in C. The full flux form operators are 
appended verbatim in E. 


3.2 General SBP Operators 

The previous fourth-order finite-domain example (as well as the sixth- and eighth- 
order examples, which are given in C) suggests that conventional SBP operators 
lead to conservative finite-domain split-form operators. To this end, we make the 
following conjecture: 

Conjecture 3.1. The discrete split-form of eq. (2.7) can be manipulated into the 
conservative form 

u t + T ,_1 Af = 0 (3.7) 

for any diagonal-norm SBP operator that can be expressed in the form of eq. (2.1) 
and for any value of the parameter a. Furthermore, the resulting local fluxes fj 
have compact support and are consistent with the original conservative flux /(it). 



The Lax-Wendroff theorem guarantees that discretely captured discontinuities are 
accurate. 

To prove this conjecture for a general banded matrix Q of halfwidth r, one would 
need to first prove that the product rule term is conservative. That is, 

(VQw + WQv) = Af e 

Because the term 'P~ 1 Qf in eq. (2.7) is already discretely conservative by eq. (A2) 
(i.e., can be manipulated into the form T’ _1 Af c ), the split spatial operator is dis- 
cretely conservative for any value of the parameter a. 


Evidence: Although a formal proof is not forthcoming at this time, considerable 
evidence suggests that the conjecture is true. Specifically, all of the banded high- 
order SBP finite-difference operators that are derived for this work (up to order 
eight) are consistent with this conjecture. 

Furthermore, consider a banded matrix Q that satisfies eq. (2.1), with elements 
[Q] (ii,i 2 ) = Q(h,i 2 ) and a halfwidth r > 1. Solving for the conservative flux f e in the 
relation (V Qw + W Qv) = Af e [25] yields the closed-form expression 

r k 

fj = ^ 1 y X w j+l v j+l-k d~ w j+l—kVj+l)q(j+l-k,.j+l) 

k=1 1=1 (3.8) 

1 <j + l, j + l-k<N 1 < j < N - 1 

fo = W lUl, J N =W N V N 

Manipulating the q^j+i-kj+i) terms produces the consistency condition 


r k 


i = EE Q(j+l-k,j+l) 

k= 1 1=1 

1 < j + l, j + l-k<N l<j<N-l 


(3.9) 


Substituting Wj = w(u) and Vj = v(u) into eq. (3.8) and using eq. (3.9) yields 


r k 

h = EE Q{j+l-k,j+ 1) wv 

k= 1 1=1 

1 <j + l,j + l-k<N 1 <j < N -1, 
Jo = wv, f N = wv 


(3.10) 


Thus, the term aQWv + (1 — a) (V Qw + W Qv) can be expressed in terms of a tele- 
scoping conservative flux, has a compact support of halfwidth r, and is consistent 
with the original conservative flux in the conservation law given by eq. (2.3). All 
of the sufficient conditions of the Lax-Wendroff theorem are met, so converged so- 
lutions using the above split operators are weak solutions to the conservation law. □ 


Remark: This conjecture has been verified symbolically for matrices with a halfwidth 
1 < r < 50. 
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3.3 A Multidomain Approach for Alternative Operators 

While great flexibility is afforded by the potential use of alternative forms of the 
equations (e.g., skew-symmetric or canonical splittings), a systematic approach is 
needed to extend alternative uniform-grid single-domain, tensor-product operators 
to complex geometries. A simple approach is to adopt a general multidomain finite- 
difference discretization that incorporates a generalized curvilinear mesh in each 
domain. 

3.3.1 Multidomain Operators 

A major obstacle in the application of high-order methods to realistic problems is 
the development a suitable grid around complex geometric features or in regions 
of strong gradients. (Constraining the grid to design-order smoothness severely 
complicates grid generation.) s equation is Multidomain techniques greatly simplify 
the grid generation process for complex configurations by breaking the geometry into 
the union of piecewise smooth quadrilateral (hexahedral) domains in two (three) 
dimensions. 

Conventional high-order SBP finite-difference techniques naturally extend to 
multidomain discretizations [21,23,26,27]. Each domain is discretized with a sta- 
ble tensor-product formulation and then connected to its adjoining neighbors using 
interface conditions that maintain the stability, accuracy, and conservation of the 
interior operators. Penalty type interface treatments, which are closely related to 
those that are used in discontinuous Galerkin and internal penalty finite-element 
approaches, are most frequently used to join the domains. Domain interfaces need 
only be Co smooth to maintain the stability, conservation and design accuracy of 
the single-domain operators. 

Existing derivations for multidomain high-order finite-difference schemes primar- 
ily focus on the divergence form of the equations. Indeed, stability and conservation 
proofs follow naturally with the use of this form. However, nothing precludes the 
use of alternative operators in each subdomain. The only additional requirement is 
the need for domain-interface coupling conditions that retain the desirable interior 
properties of the alternative operator. These interface operators must be developed 
on a case-by-case basis; however, as the examples in the next section demonstrate, 
these operators generally are simple extensions of the original coupling conditions 
introduced and extended in refs. [21,23]. 

3.3.2 Curvilinear Coordinates 

Generalized curvilinear coordinate formulations are commonly used in high-order 
formulations to facilitate the use of nonunifornr grid distributions or the definition 
of “complex” geometries. The likelihood exists that the desirable attributes of a 
Cartesian grid, alternative formulation will not survive the curvilinear mappings 
between the physical (x, y, z ) and the computational (£, 77 , Q space. However, the 
curvilinear alternative formulation can inherit the desirable Cartesian grid proper- 
ties. For example, Pirozzoli [16] recently demonstrated that pseudo-kinetic energy 


10 



can be preserved on curvilinear grids that use an alternative (skew-symmetric) for- 
mulation, (i.e., one motivated by formulations that preserve kinetic energy on a 
Cartesian grid). 


4 Examples using Burgers equation 

We use the viscous Burgers equation to validate the numerical accuracy of the alter- 
native split operators that are developed in the previous sections. In conservation 
form, the Burgers equation is 


u t + f(u) x = eu xx , f(u)=y yj , x G [0, 1], t€[0,oo) 

a o u(0,t) - eu x (0,t) - g 0 (t) = 0, aiu(l, t) - eu x (l, t) - gi(t) = 0 (4.1) 

u(0, t) + |tt(0, t)\ _ u(l,t) — \u(l,t)\ 

' “ 1 = 3 

where e represents the diffusion coefficient. Although the problem is fundamentally 
different with e = 0, we form the discrete approximation such that when e = 0 
the fully inviscid Burgers equation is recovered and only the appropriate boundary 
conditions are imposed. The boundary conditions that are given in eq. 4.1 ensure 
that the formulation is strongly well-posed; this is demonstrated further in section 
4.1. 

The Burgers equation must be split to find continuous and discrete energy esti- 
mates. In the next section, we examine numerically the difference between the full 
conservation form and the split-form by using the discretely conservative method 
that is developed herein. While our formulation allows for any value of a, only one 
a yields the energy estimate. This value of a is used in the next section to evaluate 
advantages, if any, that are gained from using the discrete formulation that satisfies 
the energy estimate. 


4.1 Energy Analysis of the Continuous Problem 

The properties and behavior of the inviscid portion of eq. (4.1) are covered in an 
energy analysis that is presented by Jameson [18]; we reproduce and supplement the 
essential results here for completeness. The energy that is given in Burgers equation 
can be analyzed directly without a linearization. If we assume a smooth solution, 
then can we split the conservation form (eq. (4.1)) by the parameter a to obtain 

U * + Q; (t) + (1 “ a ) uu x = eu xx , a; G [0,1], f G [0, oo) (4.2) 
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This equation is multiplied by u and integrated over the interval [0, L] to obtain the 
energy 


1 d_ 

2 dt 



u 


au — + (1 — a)u^u x 


dx + e J uu xx dx 
o 




dx + e / [{uu x ) x - u x u x \ dx (4.3) 


CX q 
2 ° 


+ euu x |n - / 1 - -a u u x dx - e II u x 


Equation 4.3 clearly shows that for a = 2/3 the integral term disappears and the 
energy is dependent on the boundary data and the viscous dissipation. This is re- 
ferred to as the canonical splitting or “skew-symmetric form” of the Burgers equa- 
tion [18,28]. For the viscous problem (e > 0), note that any shock that develops is 
resolvable at the continuous level and that the energy analysis remains valid. For the 
inviscid problem (e = 0), the integration bounds can be split, and the same canon- 
ical splitting parameter is found. See reference [18] for details. The well posedness 
of the boundary data is evaluated first by substituting a = 2/3 into eq. (4.3) and 
then by analyzing each boundary individually. The left boundary term (scaled by 
2 for convenience) is 



it(0, t) 


3 


2ett(0, t)u x (0, t) 


(4.4) 


Completing the squares in eq. 4.4 leads to boundary terms that can be written as 


f2 


BTn = 


\u(0,t) 


-u( 0 ,t) 2 - eu x (0,t) 


- (eu x (0,t)Y 


(4.5) 


We define 



(4.6) 


as a sensor for the hyperbolic part of the equation; do and a\ in the boundary 
conditions of eq. (4.1) serve as switches to eliminate the hyperbolic part of the 
boundary condition if waves are not propagating toward the interior of the domain 
from the boundaries. If we use the definition above, the left boundary term is 


BT 0 = — 
ao 


eu x (0, t )] 2 - [eu x { 0, f)] 2 j , a 0 = a[u(0, t)] (4.7) 


We now substitute the left boundary condition from eq. (4.1) for eu x (0,t) in eq. 
(4.7) to obtain 


BT 0 = - 
a o 


j[a 0 u(0,f) + go{t) - a 0 tt(0, t)] 2 - [d o u(0,t) - c/o(*)] 2 } • (4.8) 
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The condition ao > 0 in eq. (4.1) produces do = a, and we find that 

BT 0 = ^ - - [a 0 ix(0, t) - go(t)} 2 < (4.9) 

CLq (l Q-0 

If u{ 0, t) = a[«(0, t)] = 0, then BT$ = 0, and the energy of the continuous equation 
is bounded from above by the boundary data. If a < 0, then do = 0, which yields 
the condition 

BT 0 = — [a 0 u(0, t) + g 0 (t )] 2 - g ° ^ = ^-y [a o u(0, t) + go{t)} 2 (4.10) 

a o ao |ao| |ao| 

Again, the energy is bounded from above by the boundary data. The right boundary 
term (again scaled by 2 ) is 

2 

BT\ = —~u(L, t) 3 + 2 eu(L, t)u x (l,t) = —a\u(L, t ) 2 + 2di — 2 u(L, t)g\(t) (4.11) 

O 

where a\ = a[u(L,t)]. The same procedure is followed to show that the energy on 
the right boundary is bounded. We find that 


BTi 


^ [o-iu(L, t) + gi{t)} 2 for oi >0 

“ W\ I aiU ( L ’ ~ Slit)? for < 0 


(4.12) 


Since the energy of the continuous equation is bounded by the energy imposed 
through the boundary terms, the equation and boundary conditions in eq. (4.1) are 
strongly well-posed. 


4.2 Energy Analysis of the Semidiscrete Problem: Single Domain 

The canonical splitting can be used to construct a semi discrete operator that satisfies 
the semidiscrete analog to eq. (4.3). Using the finite-difference operators V and 
V 2 that satisfy the SBP condition on the finite domain, we discretize the skew- 
symmetric form (a = 2/3) of eq. (4.2) as 

u f = — - [DU u + UVu] + eV 2 u 

O 

+ <To'P _1 e 0 [a 0 ui - e (5u) x - g 0 (t)] 

+ a 1 V~ 1 e 1 [aiUAr - e (<Su ) N - 51 (f)] (4.13) 

U = diag(u), e 0 = (1,0,...) T , ei = (..., 0, 1 ) T 
rti + kil _ un - [un\ 
a ° ~ 3 ai “ 3 

The proper values of 00 and a \ are determined in the following energy analysis. Both 
forms of V 2 in eq. (Dl) (see appendix D) satisfy the energy analysis, but we use 
the narrow stencil here for illustrative purposes. The simultaneous approximation 
term (SAT) penalty method [29] is used in eq. (4.13) to satisfy the mixed boundary 
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conditions with specified boundary data. For the finite domain, a diagonal norm 
V is chosen to ensure that matrix multiplication of the norm and the nonlinear 
coefficient matrix commutes. This choice reduces the overall accuracy for hyperbolic 
problems to (p + 1), where 2 p is the internal accuracy and p is the accuracy at the 
boundary. For parabolic problems, the narrow- stencil second derivative operator 
V 2 = V-\n + BS) constructed with the diagonal norm, provides a global (p + 2) 
order of accuracy [30]. Premultiplying eq. (4.13) by u T V and using eq. (2.1) yields 
the discrete energy 


x“T _ ll u llp = — ^u T (QU + UQ) u + €u t TZvl + eu T £>Su 
+ u i e 0 {o- 0 [a 0 ui - e(5u) 1 -£f 0 (i)]} 

+ u 1 ei {<7i [aiu N - e (<Su) N - gi(t)}} 

Adding eq. (4.14) to its transpose yields 
d 2 

— ||u|||> = — -u t BUu + 2eu T lZu + 2eu T BSu 

+ 2tticr 0 ( a 0 ui - e [<Su) 1 - g 0 (t)] 

+ 2u n g\ {d\u N - e[Su) N - gi (t)] 

Simplifying with eq. (2.1) yields 


(4.14) 


(4.15) 


U 


dt 


v 


= 2err TZu 


+ 2u\ \ <7 0 [a 0 ui - 5o(i)] + 2 U 1 


+ 2u n { <7i [fil'd tv - gi(t)j - -u 2 n 


2eu\ [(<j 0 + 1) (<Su) 1 ] 

+ 2eu N [(-<7i + 1) {Su) N j 


(4.16) 


Above, 1Z is negative semi definite and ensures that the energy only decays. The 
viscous boundary terms cancel for <7o = — 1 and oq = 1, which leaves 

— ||u||?> = 2eu 7 IZu + -BTo + BT\ 
dt 

BT 0 = ui [a(ui)ui - 2d 0 ui + 2g 0 (t)} ( 4 - 17 ) 

BT i = utv [2fiidjv — a(ujv)'djv — 2g\(t)\ 


For the discrete equation to be strongly well-posed, the discrete energy must be 
bounded from above by the imposed boundary data. For the left boundary when 
a(tti) > 0 and using completing the squares, 

—— a 2 1 

BT 0 = —au\ + 2uig 0 = — (aui - g 0 ) 2 . (4.18) 

a a 

The result in eq. 4.18 is equivalent to continuous result found in eq. (4.9). 
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Completing the squares at the left boundary if a(ui) < 0 yields 

... 2 2 

BT 0 = au\ + 2u 1 g 0 = - — {aui + g 0 )' 2 (4.19) 

|a| |a| 

Again, note the similarity between eq. 4.19 and the continuous result given in eq. 
(4.10). The right boundary also mimics the continuous case, so the semidiscrete 
equation is strongly well-posed. Note that the same penalty is used for the inviscid 
and the viscous Burgers equations and that the equation is well-posedness in both 
cases. 


4.3 Energy Analysis of the Semidiscrete Problem: Multidomain 

Often multiple semidiscrete domains can be used advantageously to discretize a 
conservation law. Interface penalties are derived such that the semidiscrete energy 
is unaffected by the exchange of information between two domains. For a two domain 
problem on the interval x £ [—1,1] divided at x = b. we use the discretization 


x« = 


/ y * ' 7 ry » ' 7 / y * ' 7 I 

> x 2 :■■■■> X 7V 2 -1’ X TV 2 ’ 


4 1} = -1 + (z - 1) » * = l,2,...,Ni 

x? = b + (i - 1) ^ , i = l,2,...,N 2 


where N± and N 2 are the number of uniform cells in each of the two domains. The 
solution on each domain is denoted by and u^ 2 ), respectively. Operators that are 
defined on each domain are denoted similarly. The semi discrete Burgers equation 
for the multidomain problem is 


u| 1] = (£> (1) t/ (1) u (1) + t/ (1) P (1) u (1) ) +el4 1) u (1 > 

+ (P (1) ) | - koiuf^ + d 0 i (s (1) u (1) )^ - (5 (2 ) u (2) ^ | 

+ (V (i) ) (<S (i) ) e i 1] ^io (^jv! - «i 2) ) 

u{ 2) = ~ (p>( 2 )[/( 2 ) u (2) + t/ (2) P (2) u (2) ) + eP$ 2) u (2) 

+ (V (2) ) e[ ) 2 ' > - k w u ^*) +er 0 i (^5 (2, u (2) ^ - (s (1) u (1) )^ | 

+ (V (2) ) (cS (2) ) e^ 2) eno (*4 2) - “tv!) 

C/ (1) = diag(u (i) ), C/ (2) = diag(u (2) ), = (..., 0, 1) T , e^ 2) = (1, 0, . . . ) T 

(4.20) 

where the boundary condition penalties have been dropped for clarity of the inter- 
face treatment. The energy is calculated by premultiplying the two equations by 
(u (!)) pW anc l ( u ( 2 )) p( 2 ) respectively. The equations are then added to their 
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transpose as 



(4.21) 


The terms that are related to the boundary data have been removed. Because we 
want to construct the interface penalty such that it conserves the energy, we must 
choose the free parameters such that the right-hand side of eq. (4.21) vanishes. This 
results in the following constraints: 


koo = 2 U Nl > 
ho = —loi — T 


ho = -koi, k n = -^u { {‘ ) 
m = ^oi + l> no = —^oi 


The formal order of accuracy of the penalty terms must be the same as the order of 
the boundary finite-difference approximations. This requires that the value of &oi 
be a function of u. We choose an equal weighting of the values on either side of the 
interface: 


*01 


u 


( 1 ) 

Ni 


+ «< 2) 


(4.22) 


For the viscous terms, foi = yields an equal weighting of all of the viscous 
penalty terms; thus, this is the chosen value. The full semidiscrete set of equations 
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for the two-domain problem is 


1 

u, = - 


.(!) _ 


u (1 )) + eP^u (1 ) 
(^cS (1, u (i) j i - g 0 (t) 


+ 


-1 


.(!) 


.(!) 




f . „(2) ' 

( 1 ) ( 1 ) U N 1 + U 1 ( 2 ) 
U N 1 U N 1 9 U 1 


(s« u«) w -(s^V 2 ^ 

-1 / ,.,\T 


(pW) (S (1) ) ej 1} ( 


(i) .,(2) 


u n[~ u 


i 1 ’ 


uf } = - * (p( 2 )[/(2) u (2) + C/ (2) P (2 V 2) ) + eP^ 2) u (2 ) 


-i 


,( 2 ) 


? , (1 ) , .,( 2 ) 

(2) (2) W 7Vi + U 1 (1) 

u\ 'u\ — z v> N [ 


n “i 


(4.23) 


€ 

~2 

e 

+ 2 


+ 


( 5 ( 2 ) u ( 2 )) -(^uW 

(V (2) ) _1 (^ (2) )^ 


Ah 


(2) / ,(2) 


'«! — U 


-1 


,( 2 ) 


~ ( 2 ) ( 2 ) 

°1 U 7Va - 6 




In this form, the boundary penalty does not contribute to the energy. The en- 
ergy changes, however, as a result of the one-sided stencils at the interface. The 
semidiscrete equations remain strongly well-posed. 


5 Numerical Tests 


The new discrete operators were implemented for the Burgers equation in a one- 
dimensional finite-difference solver with a uniform grid. Integration in time was 
conducted with a five-step, fourth-order Runge-Kutta scheme [31]. The step size in 
time was chosen such that the temporal error is negligible in comparison with the 
spatial error. 

Two problems are evaluated. For the first problem, the accuracy of the three 
operator sets that are developed herein is tested by using the skew-symmetric ( a = 
2/3) and the conservation (a = 1) forms of the viscous Burgers equation. The 
second problem tests the accuracy of the multidomain interface closures that are 
developed based on the new operator sets. 


17 



5.1 Accuracy Validation 

The viscous Burgers equation is used to test the accuracy of the new operator when 
it is applied to the nonlinear equation. The test problem is defined as 

u f = V~ l A (af c + (1 - a)f e + fy) 

+ ooT^e o [—a,QU\ - e (5u) 1 - g 0 (t)] 

+ ai'P _1 ei [ Qj\u N - e (Su ) N - gi(t)] 

9 o{t ) = — + 1 ' t ^ v(-l,t) - ev x (-l,t) (5.1) 

9i(t) = ? - ev x (l,t) 

A.x 

v(x,t) = 1 , xG[xi,xjv], te[0,T] 

x ' + 2t + w 

where e = 1 and the initial condition is uo = u(x,0). The definition of the viscous 
flux fy is given in D. The values for the penalty coefficients a are those that are 
identified in section 4. Note that the sign of the convective term has changed for 
this problem. This is to satisfy the exact solution v(x,t). The problem is simulated 
on the interval [—1, 1] up to T = 0.05. Five grids are used to determine the order 
of accuracy of the set of operators. The L 2 and L°° error norms of the conservation 
form and the canonical split-form are compared in tables 1 through 3. Figure 1 
shows the numerical solution for N = 64 at different times. All simulations use the 
narrow-stencil viscous operator. 


Table 1. Error Norms and Convergence Rates for Viscous Burgers Equation Using 
Canonical Split and Conservation Forms, and Conservative Narrow-stencil (2-4-2) 


Viscous Operator. 



Canonical Split 

Conservation 

N 

L' z Error 

L z Rate 

L°° Error 

L°° Rate 

L z Error 

L z Rate 

L°° Error 

L°° Rate 

32 

1.64e-03 

- 

2.70e-03 

- 

1.29e-03 

- 

2.44e-03 

- 

64 

1.10e-04 

3.90 

1.91e-04 

3.82 

8.86e-05 

3.86 

1.75e-04 

3.80 

128 

6.97e-06 

3.98 

1.21e-05 

3.97 

5.66e-06 

3.97 

1.17e-05 

3.90 

256 

4.38e-07 

3.99 

7.66e-07 

3.99 

3.56e-07 

3.99 

7.38e-07 

3.99 

512 

2.74e-08 

4.00 

4.79e-08 

4.00 

2.23e-08 

4.00 

4.62e-08 

4.00 


Table 2. Error Norms and Convergence Rates for Viscous Burgers Equation Using 
Canonical Split and Conservation Forms, and Conservative Narrow-stencil (3-6-3) 


Viscous Operator. 



Canonical Split 

Conservation 

N 

L 2 Error 

L z Rate 

L°° Error 

L°° Rate 

L 2 Error 

L 2 Rate 

L°° Error 

L°° Rate 

32 

2.34e-04 

- 

4.62e-04 

- 

2.83e-04 

- 

6.24e-04 

- 

64 

4.56e-06 

5.69 

9.15e-06 

5.66 

5.61e-06 

5.65 

1.30e-05 

5.58 

128 

7.52e-08 

5.92 

1.51e-07 

5.92 

9.31e-08 

5.91 

2.19e-07 

5.90 

256 

1.19e-09 

5.98 

2.40e-09 

5.98 

1.48e-09 

5.98 

3.49e-09 

5.97 

512 

2.00e-ll 

5.90 

3.98e-ll 

5.92 

2.33e-ll 

5.99 

5.37e-ll 

6.02 


The expected order of accuracy is p + 2, where p is the boundary accuracy. 
This design order is achieved for both the skew-symmetric and conservation forms 
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Table 3. Error Norms and Convergence Rates for Viscous Burgers Equation Using 
Canonical Split and Conservation Forms, and Conservative Narrow-stencil (4-8-4) 


Viscous Operator. 



Canonical Split 

Conservation 

N 

L' z Error 

L 2, Rate 

L°° Error 

L°° Rate 

L z Error 

L z Rate 

L°° Error 

L°° Rate 

32 

1.52e-04 

- 

3.64e-04 

- 

1.89e-04 

- 

4.42e-04 

- 

64 

6.97e-07 

7.77 

2.17e-06 

7.39 

8.58e-07 

7.78 

2.46e-06 

7.49 

128 

3.85e-09 

7.50 

1.78e-08 

6.93 

4.45e-09 

7.59 

1.92e-08 

7.00 

256 

2.74e-ll 

7.13 

1.91e-10 

6.55 

2.94e-ll 

7.24 

1.98e-10 

6.60 

512 

1.74e-12 

3.98 

2.65e-12 

6.17 

1.73e-12 

4.09 

2.69e-12 

6.20 



Figure 1. Numerical approximation to skew-symmetric form of viscous Burgers 
equation with 64 uniform cells and the (4-8-4) narrow stencil conservative operators. 


for the (2-4-2) operator. However, the (3-6-3) and (4-8-4) operators exhibit greater 
accuracy than design order convergence. The (2-4-2) operator does not deviate from 
design order because the maximum order that can be achieved is equivalent to the 
design order. This is not the case for the (3-6-3) and (4-8-4) operators. For these 
operators, the accuracy is at least design order for well resolved solutions. For the 
(2-4-2) operator, the conservation form is slightly more accurate than the canonical 
split-form. For the (3-6-3) and (4-8-4) operators, the canonical split-form is more 
accurate. 
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5.2 Multidomain Validation 


The test problem that was used to test the interfaces is described by the exact 
solution 

v(x, t) = 1 — tanh , x G [—10, 10], t G [0, T] (5.2) 

where e = 0.25, xq = —5.0, and the initial condition is uo = v(x, 0). The problem is 
simulated up to T = 10.0 with four grid resolutions. The skew-symmetric form of 
the discrete equations is used. 

The L 2 and L°° errors are tabulated in tables 4 through 6 for each operator on 
(1) a single domain, (2) two domains with equivalent spacing, and (3) two domains 
with grid spacing that changes discontinuously at the interface by a factor of r = 2. 
(Herein, the jump in spacing is denoted as the “compression ratio.”) The interface 
between domains is located at x = 0.0. The discrete solution with 256 uniform cells 
on a single domain for the (2-4-2) operator set is shown in figure 2. 



Figure 2. Discrete solution to viscous Burgers equation using 256 equispaced nodes 
with (4-8-4) narrow operators. 

Tables 4 through 6 demonstrate that the interface has a negligible effect on the 
accuracy for this well-resolved test case. In addition, the treatment has been verified 
to be design order accurate, even when a discontinuous grid resolution is used across 
the subdomains. 
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Table 4. Error Norms and Convergence Rates of (2-4-2) Operator for Single Domain, 
Two Domains with Uniform Resolution, and Two Domains with compression ratio 


r = 2. 


Single Domain 


Resolution 

Error 

N 

L z Error 

L z Rate 

L°° Error 

L°° Rate 

64 

8.89e-03 

- 

1.13e-02 

- 

128 

6.37e-04 

3.80 

8.86e-04 

3.68 

256 

4.19e-05 

3.93 

6.17e-05 

3.84 

512 

2.66e-06 

3.98 

4.04e-06 

3.93 


Two Domains, r = 1 


Resolution 

Error 

iVi 

n 2 

L 2 Error 

L z Rate 

L°° Error 

L°° Rate 

32 

32 

8.91e-03 

- 

1.14e-02 

- 

64 

64 

6.38e-04 

3.80 

8.87e-04 

3.68 

128 

128 

4.19e-05 

3.93 

6.17e-05 

3.85 

256 

256 

2.66e-06 

3.98 

4.04e-06 

3.93 


Two Domains, r = 2 


Resolution 

Error 

N\ 

n 2 

L 2 Error 

L z Rate 

L°° Error 

L°° Rate 

32 

64 

6.26e-04 

- 

8.65e-04 

- 

64 

128 

4.12e-05 

3.93 

6.08e-05 

3.83 

128 

256 

2.61e-06 

3.98 

3.97e-06 

3.94 

256 

512 

1.64e-07 

4.00 

2.50e-07 

3.99 


Table 5. Error Norms and Convergence Rates of (3-6-3) Operator for Single Domain, 
Two Domains with Uniform Resolution, and Two Domains with compression ratio 


r = 2. 


Single Domain 


Resolution 

Error 

N 

L z Error 

L z Rate 

L°° Error 

L°° Rate 

64 

3.48e-03 

- 3.96e-03 

- 


128 

9.88e-05 

5.14 

1.33e-04 

4.89 

256 

1.89e-06 

5.71 

3.23e-06 

5.37 

512 

3.13e-08 

5.92 

5.41e-08 

5.90 


Two Domains, r 

= 1 


Resolution 

Error 

Ni N 2 

L z Error 

L z Rate 

L°° Error 

L°° Rate 

32 32 

3.47e-03 

- 

3.91e-03 

- 

64 64 

9.87e-05 

5.14 

1.33e-04 

4.88 

128 128 

1.89e-06 

5.71 

3.21e-06 

5.37 

256 256 

3.11e-08 

5.92 

5.37e-08 

5.90 


Two Domains, r 

= 2 



Resolution 

Error 

TVi 

n 2 

L z Error 

L z Rate 

L°° Error 

L°° Rate 

32 

64 

9.82e-05 

- 

1.34e-04 

- 

64 

128 

1.87e-06 

5.72 

3.15e-06 

5.41 

128 

256 

3.11e-08 

5.91 

5.05e-08 

5.96 

256 

512 

5.54e-10 

5.81 

7.55e-10 

6.07 


6 Conclusion 

We have developed a class of split-form finite- difference operators that satisfy the 
sufficient conditions of the Lax-Wendroff theorem. These operators are applicable to 
the conservation or divergence form of the conservation law, but facilitate operator 
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Table 6. Error Norms and Convergence Rates of (4-8-4) Operator for Single Domain, 
Two Domains with Uniform Resolution, and Two Domains with compression ratio 


r = 2. 


Single Domain 


Resolution 

Error 

N 

L z Error 

L z Rate 

L°° Error 

L°° Rate 

64 

1.97e-03 

- 

2.00e-03 

- 

128 

2.88e-05 

6.10 

3.56e-05 

5.81 

256 

1.81e-07 

7.31 

3.37e-07 

6.72 

512 

8.53e-10 

7.73 

1.36e-09 

7.95 


Two Domains, r = 1 


Resolution 

Error 

iVi 

n 2 

L z Error 

L z Rate 

L°° Error 

L°° Rate 

32 

32 

2.16e-03 

- 

2.70e-03 

- 

64 

64 

2.89e-05 

6.23 

3.79e-05 

6.15 

128 

128 

1.84e-07 

7.29 

2.93e-07 

7.02 

256 

256 

1.01e-09 

7.51 

1.67e-09 

7.46 



Two Domains, 

r = 2 



Resolution 

Error 

iVi 

n 2 

L z Error 

L z Rate 

L°° Error 

L°° Rate 

32 

64 

3.01e-05 

- 

4.44e-05 

- 

64 

128 

2.67e-07 

6.82 

4.36e-07 

6.67 

128 

256 

3.43e-09 

6.28 

4.63e-09 

6.56 

256 

512 

1.09e-10 

4.98 

4.83e-ll 

6.58 


splitting at the discrete level to improve accuracy or robustness, depending on the 
conservation law. The method for constructing a conservative split operator is 
illustrated for the (2-4-2) operator with a general splitting parameter. Higher order 
operators were also derived and are supplied in an accompanying text file. We also 
developed fully conservative operators for linear viscous terms. 

The specific operators that are derived herein were tested on the conservation and 
canonical split forms of Burgers equation. Energy stability of the discrete form was 
proven, including boundary and interface closures. The split and conservation forms 
yield very similar results for Burgers equation. The solutions converged at least at 
the design order of accuracy in all cases, even when multiple domains were used 
with discontinuous grid resolutions across interfaces. The interface treatment that 
was tested herein had a negligible effect on the accuracy of the solution. Nonlinear 
dissipation operators for split forms will be presented in a future paper. 
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Appendix A 


Discrete Conservation 

The summation-by-parts operators are discretely conservative, as demonstrated 
in the following lemma. 

Lemma A.l. All differentiation matrices that satisfy the SBP convention given in 
eg. (2.1) are discretely conservative in the V-norm. 

Proof. Consistency of V implies that all rows in V the matrix sum to zero. Because 
V is nonsingular, the following matrix relations hold: 

V sb pl = V- 1 Ql = Q1 = 0 , 1 = (1, ■ ■ ■ ,1) T , 0 = (0, • • • , 0) T (Al) 

which demonstrates that all rows of Q sum to zero. Forming the dotproduct of 
Q + Q T with the unit vector and combining the result with that obtained in eq. 
(Al) yields 

(Q + Q T ) 1 = Q T 1 = B1 = (-1,0, . . . , 0, 1) T (A2) 

Taking the transpose of eq. (A2) and combining it with the definition of an SBP 
operator yields 

1 t Q= l T VV sbp = (-1,0,..., 0,1) (A3) 

Thus, eq. (A3) demonstrates that the columns of T> s b p are conservative in the V 
norm (i.e., all interior columns of the matrix V sum to zero, and the left and right 
boundary columns sum to —1 and 1, respectively, after premultiplication by V) . □ 
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Appendix B 


Split Continuous Equations 


To utilize a split operator, the flux / in eq. (2.3) must be the product of at least 
two functions: f(u) = v(u)w(u). Equation (2.3) is split by using the product rule 

u t + af(u) x + (1 - a) [v(u)w(u) x + w(u)v(u ) x ] = 0. (Bl) 

This splitting does not change the weak solution to the governing equation, nor does 
it alter the Rankine Hugoniot relation. To show this, we start with the conservation 
law in one dimension: 

xr 

d f 

— udx+ f(u)\l* = 0, xe[x L ,x R }, te[ 0,oo) (B2) 

X L 

with a discontinuity located at X&. At this point, note that neither eq. (2.3) nor 
eq. (Bl) is valid across the discontinuity. Both differential forms are valid on either 
side of the discontinuity, with the mild restriction that w(u) x and v{u) x exist in all 
regions where f(u) x exists. Equation (B2) can be rewritten [32] as 


d_ 

dt 


I udx+ 


XR 

J t J udx+ f{u ) \ x x r l 


= 0 


XL 


X 


+ 

d 


where u is smooth over the interval of integration. If we apply the Leibniz integral 
rule, then 


d 

s 


XL 


ut dx + Su(x d 


XR 

J' 


ut dx — Su(x d 


+ /(«)!£ = o 


where S = ^ is the propagation speed of the discontinuity. Either eq. (2.3) or eq. 
(Bl) can be used to substitute for u t , because the two forms are equal for smooth 
solutions. The conservation law becomes 


X d XR 

J —f{u) x dx + Su(x d ) + J -f(u) x dx - Su(xj) + f(u) \ x x * = 0 

XL 4 

Upon simplification, we get the Rankine Hugoniot relation 

[f( x d ) - f( x d )] - 5 H x d ) - u ( x d )] = 0 ( B3 ) 

Thus, both eq. (2.3) and eq. (Bl) satisfy the sufficient condition of the Lax-Wendroff 
theorem which specifies that the continuous weak solution of the discretized equa- 
tions must satisfy the Rankine Hugoniot relations. 
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Appendix C 


High-Order Skew-Symmetric Fluxes 


C.l The Fourth-Order Operator: (3-6-3) 

The full (3-6-3) operator is available upon request. The use of electronic format ne- 
cessitates the brevity of this presentation and easier implementation. For reference, 
the interior flux operator is 



/ Ui-2 

\ 

T 

( 2 

0 

0 

1 

0 

0 ^ 


( Ui-2 \ 


Ui - 1 



0 

-16 

0 

-9 

1 

0 


Ui - 1 

1 

Ui 



0 

0 

74 

45 

-9 

1 


Ui 

360 

^2+1 



1 

-9 

45 

74 

0 

0 


^2+1 


^2+ 2 



0 

1 

-9 

0 

-16 

0 


^2+2 


\ u i+3 

/ 


V o 

0 

1 

0 

0 

2 ) 


\ u i + 3 / 


(Cl) 


C.2 The Complete Fifth-Order Operator: (4-8-4) 

The full operator can be found in the accompanying text file. The flux for the 
interior is 


ft = 

J 840 


Ui - 3 

\ 

T 

f 

-l 

0 

0 

0 

1 

2 

Ui -2 




0 

29 

3 

0 

0 

Iff 

3 

Ui - 1 




0 

0 

139 

3 

0 

-28 

Ui 




0 

0 

0 

533 

3 

112 

Ui +1 




1 

2 

16 

3, 

-28 

112 

533 

3 

Ui +2 




0 

1 

2 

16 

3, 

-28 

0 

Ui + 3 




0 

0 

1 

2 

16 

3, 

0 

Ui + 4 

J 


{ 

0 

0 

0 

1 

2 

0 


_ i 

if 

3 

-28 

0 

139 

3 

0 

0 


O 

O 


l Ui - 3 \ 

0 0 


Ui -2 

o ^ 
1 - 


Ui - 1 

3 2 


Ui 

0 0 


Ui+1 

0 0 


u i +2 

f o 


Ui + 3 

0 -1 J 


\ U i+ 4 J 


( C 2) 
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Appendix D 


Conservative Finite-Difference Viscous Terms 


Many conservation laws of interest also have viscous terms. We show that for the 
linear viscous terms eu xx the summation-by-parts narrow-stencil approximation to 
the second derivative can be cast into a conservative form. The SBP finite-difference 
operators for the second derivative can be written as 

V -2 = V 2 or V -2 = V _1 (1Z + BS), 1Z = 1Z T , v T lZv < 0 (Dl) 

where S represents a first derivative operator. In eq. ??, T > 2 = T> 2 is called the wide- 
stencil approximation to the second derivative, which requires that S = V; V 2 = 
V _1 ( TZ + BS ) is the narrow-stencil approximation to the second derivative and does 
not impose the same restriction on S [30]. The narrow-stencil approximation has 
superior accuracy and stability properties and should be used where possible. 

The conservation form of the narrow stencil is found by constructing an approx- 
imation to the viscous flux, eu x , at the flux points such that 

M v = -eA^u = - e (K + BS) u (D2) 

This is solved by using a linear system for an appropriate TZ and BS that satisfies 
the SBP condition. The viscous flux 


fy = —eTZu 


(D3) 


yields a fully conservative discrete approximation to —eT> 2 U = Af y. The operator 
TZ for the 2p = 4 scheme is presented herein, while the TZ operators for 2p = 6 and 
2p = 8 are appended verbatim in E. For the fourth-order case, we begin with the 
existing fourth-order operator of Mattsson et al. [30], which can be expressed with 
the relation V 2 W = A7?,u = (7 Z + BS) u, where 


TZ5x 


5# 

4 §9 


48 

0 

59 

yj 

0 

1 

48 1 

5§ 4 

4 §5 
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48 
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4 §9 

P 

48 

0 

48 1 

i 4 

% 

12 
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2 
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0 

0 



0 

0 

0 

0 

1 

12 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

j_ 

P 
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0 

0 
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0 

0 

0 

0 
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\ 

0 

0 

0 

0 
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(D4) 
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and 


BSSx = 


( ¥ 

0 


0 

V o 


-3 


3 
2 

0 0 


0 0 0 
0 0 0 0 


0 0 
0 0 


0 0 
0 0 


0 0 0 0 
0 0 0 0 


0 0 0 

1 -3 


Solving for 1Z in the relation AIZu = (1Z + BS ) u yields 

/ 
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(D5) 


(D6) 


Note that this operational form is only design order accurate for a constant diffu- 
sivity e. 

A similar manipulation can be performed by using the sixth- and eighth-order 
diffusion operators. These operators are appended verbatim in E. 
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Appendix E 


Stencil Coefficients 

The following section includes the coefficients for implementation of the split 
form conservative finite difference simulations. The given coefficients can be inserted 
directly into source code using “cut-and-paste” . 

E.l Definitions 

• First derivative of vector u, D\u = V~ 1 Qu, [see eq. (2.9)] . 

• Second derivative of vector u, D^u = + BS)u, [see eq. (3.1)]. 

• Split first derivative of w(u)v(u), 

aDiWv + (1 - a){WDiv + VD lW ) = T~ l A f 
, [see eq. (2.24)]. 

• Conservative Second Derivative, V~ l A f = —eV~ 1 A'Ru, [see eq. (3.2)]. 

For each scheme [ (2-4-2), (3-6-3), and (4-8-4) ], the following are given: 

• The flux points for the local stencil, written in terms of a general flux / = v*w. 

• The diagonal of the P-norrn V scaled by the grid spacing 

• The diagonal of the inverse of the P-norrn V scaled by the inverse of grid 
spacing 

• The nonzero entries of the Q matrix for the SBP first derivative operator (for 
reference) 

• The nonzero entries of the matrix 7Z, used to calculate the conservative flux 
terms 

• The nonzero entries of the matrix 7Z, which approximates the second derivative 

• The nonzero entries of the matrix B * S, which yields the boundary differen- 
tiation operator 

Note that only the first s+1 entries are shown. The (s+l)th entry is the interior 
operator, and just shifts indices for subsequent entries in the matrix. The respective 
symmetry or skew symmetry properties of each matrix is as follows: 

• The flux points are per-symmetric. Replace Vi and Wj with ujv-(i-i) and 
wjv-D-i)) respectively, at right boundary i=l,2,...,N. Note that for Burgers 
equation im = u t and Vj = Uj/2. 
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• The P-norm V is per-symmetric. 


Pi,i Pn~ ( i— 1),N— (i— 1) j * 1) 2, N 


The inverse P-norm V 1 is per-symmetric. 


■p - 1 — -p^ 1 j — io /V 


The boundary block portion of Q is per-skew-synrnretric. 

QN-(i-l),N-(j-l) = = 1) 2, s, j = 1,2, s 

The boundary block portion of 7Z is per-skew-synrnretric. 

P'N—(i),N—(j—l) P"i)j ’ 4 0, 2, ..., s — 1, j 1,2, ..., s 

The boundary block portion of 1Z is per-symmetric. 

P-N— (i— 1),N— (j— 1) P’i,j i 4 1, 2, ..., S, J 1,2, ..., S 


• The boundary contribution 1? * 5 is per-symmetric. 

BSN,N-(j-i) = & * S\ j,j = 1, 2, ..., A 4 


E.2 Coefficients 

E.2.1 (2-4-2) Scheme Skew- Symmetric 

Flux point j: 

= -0.08333333333333333333333333333*a+0.0000000000000000000000000000 
Vj-i*w j+0 = 0.0000000000000000000000000000*a;+0.0000000000000000000000000000 
Vj-^Wj+i = 0.08333333333333333333333333333*a— 0.08333333333333333333333333333 
Vj-i*w j+2 = 0.0000000000000000000000000000*o+0.0000000000000000000000000000 
v j+0 * Wj -i = 0.0000000000000000000000000000*o+0.0000000000000000000000000000 
v j+0 *Wj+o = 0.5833333333333333333333333333*0+0. 0000000000000000000000000000 
Vj+o*Wj+\ = —0.6666666666666666666666666667*a+0.6666666666666666666666666667 
v j+0 *w j+2 = 0.08333333333333333333333333333*o-0.08333333333333333333333333333 
++ i*''+‘-i = 0.08333333333333333333333333333*o-0.08333333333333333333333333333 
Vj+i*Wj + o = —0.6666666666666666666666666667*0+0.6666666666666666666666666667 
v j+1 *w j+ i = 0.5833333333333333333333333333*0+0. 0000000000000000000000000000 
v j+1 *w j+2 = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 

v j+ 2 *w.j-i = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
Vj + 2 *w j+ o = 0.08333333333333333333333333333*a— 0.08333333333333333333333333333 
v j+ 2 *w j+ i = 0.0000000000000000000000000000*o+0.0000000000000000000000000000 
v j+2 *w j+2 = —0.08333333333333333333333333333*0+0. 0000000000000000000000000000 
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Flux point j=3 

vi*wi = -0.03125000000000000000000000000*a+0.0000000000000000000000000000 
V! *w 2 = 0.0000000000000000000000000000* a + 0.0000000000000000000000000000 
Vi *w 3 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
V1*W4 = 0.03125000000000000000000000000*a— 0.03125000000000000000000000000 
V! *w 5 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
v 2 *wi = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 2 *w 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 2 * W3 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 2 * U , A = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 2 * w 5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 3 * Wl = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 3 *w 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 3 * xij 3 = 0.5312500000000000000000000000 * a + 0.0000000000000000000000000000 
v 3 *w A = —0.6145833333333333333333333333*a+0. 6145833333333333333333333333 
v 3 *u> 5 = 0.08333333333333333333333333333*a— 0.08333333333333333333333333333 
v A *w 1 = 0.03125000000000000000000000000*a— 0.03125000000000000000000000000 
v A * W 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v A *w 3 = —0.6145833333333333333333333333*a+0. 6145833333333333333333333333 
VA * WA = 0.5833333333333333333333333333 * a + 0.0000000000000000000000000000 

V A * W 5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

V 5 * Wl = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 5 *w 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 5 *w 3 = 0.08333333333333333333333333333*a— 0.08333333333333333333333333333 
v 5 * xu A = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 5 *w 5 = —0.08333333333333333333333333333*a+0. 0000000000000000000000000000 


Flux point j=2 

vi*W! = —0.1145833333333333333333333333*a+0. 0000000000000000000000000000 
vi *w 2 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
Vl *w 3 = 0.08333333333333333333333333333*a— 0.08333333333333333333333333333 
vi*W4 = 0.03125000000000000000000000000*a— 0.03125000000000000000000000000 
v 2 *w 1 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 2 *w 2 = 0.6145833333333333333333333333 * a + 0.0000000000000000000000000000 
v 2 *w 3 = —0.6145833333333333333333333333*a+0. 6145833333333333333333333333 
xj 2 * w A = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 3 *w { = 0.08333333333333333333333333333*a— 0.08333333333333333333333333333 
v 3 *w 2 = —0.6145833333333333333333333333*a+0. 6145833333333333333333333333 
V3 * xv 3 = 0.5312500000000000000000000000 * a + 0.0000000000000000000000000000 

XJ 3 * w A = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

v A *w i = 0.03125000000000000000000000000*a— 0.03125000000000000000000000000 
V 4 *w 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
VA * w 3 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v A *w A = -0.03125000000000000000000000000*a+0. 0000000000000000000000000000 
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Flux point j=l 

vi*wi = 0.5000000000000000000000000000* a + 0.0000000000000000000000000000 
v 1 *w 2 = —0.6145833333333333333333333333*a+0. 6145833333333333333333333333 
V!*w 3 = 0. 08333333333333333333333333333*a— 0.08333333333333333333333333333 
vi*W4 = 0.03125000000000000000000000000*o— 0.03125000000000000000000000000 
v 2 * Wl = —0.6145833333333333333333333333*a+0. 6145833333333333333333333333 
v 2 *w 2 = 0.6145833333333333333333333333 * a + 0.0000000000000000000000000000 
v 2 *w 3 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

V 2 * w 4 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

v 3 *w i = 0.08333333333333333333333333333*a— 0.08333333333333333333333333333 
v 3 * w 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 3 *w 3 = —0.08333333333333333333333333333*a+0. 0000000000000000000000000000 
v 3 * WA = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v A * Wl = 0.03125000000000000000000000000*a— 0.03125000000000000000000000000 
v A *w 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 4 *w 3 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v A *w A = -0.03125000000000000000000000000*a+0. 0000000000000000000000000000 


Flux point j=0 

Vi*wi = 0.0000000000000000000000000000 * a + 1.000000000000000000000000000 


The P-norm V is : 

Pi,i/dx = 0.35416666666666666666666666666666666667 
p 2j2 /dx = 1.2291666666666666666666666666666666667 
p 3 , 3 /dx = 0.89583333333333333333333333333333333333 
p A , A /dx = 1.0208333333333333333333333333333333333 
Pb,b/dx = 1.0000000000000000000000000000000000000 

The inverse P-norm V~ l is: 

■P M ' * dx = 2.8235294117647058823529411764705882353 
V 2 2 * dx = 0.81355932203389830508474576271186440678 
V 3 3 * dx = 1.1162790697674418604651162790697674419 
V A \ * dx = 0.97959183673469387755102040816326530612 
* dx = 1.0000000000000000000000000000000000000 


The Qij * dx nonzero terms are: 

Qi,i = -0.50000000000000000000000000000000000000 
Qi’ 2 = 0.61458333333333333333333333333333333333 
Qi’ 3 = -0.083333333333333333333333333333333333333 
Q 1A = -0.031250000000000000000000000000000000000 
Q 2l = -0.61458333333333333333333333333333333333 
22^3 = 0.61458333333333333333333333333333333333 
Q 3 .i = 0.083333333333333333333333333333333333333 
Q 3 ’ 2 = -0.61458333333333333333333333333333333333 
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03 . 4 = 0.61458333333333333333333333333333333333 
Q 3 ’ 5 = -0.083333333333333333333333333333333333333 
Ot’i = 0.031250000000000000000000000000000000000 
g 4 ’ 3 = -0.61458333333333333333333333333333333333 

24.5 = 0.66666666666666666666666666666666666667 
g 4 ’ 6 = -0.083333333333333333333333333333333333333 
Q 5 ’ 3 = 0.083333333333333333333333333333333333333 

25.4 = —0.66666666666666666666666666666666666667 

25.6 = 0.66666666666666666666666666666666666667 
2 5 ’ 7 = -0.083333333333333333333333333333333333333 

The (2-4-2) viscous terms. 

The TZi -j * dx nonzero terms: 

72 0 .1 = -1.8333333333333333333333333333333333333 
n 02 = 3.0000000000000000000000000000000000000 
TZ 0:3 = -1.5000000000000000000000000000000000000 
n 0 A = 0.33333333333333333333333333333333333333 
72 u = -1.1250000000000000000000000000000000000 
n h2 = 1.2291666666666666666666666666666666667 
72i’ 3 = -0.083333333333333333333333333333333333333 
72 - i ’ 4 = -0.020833333333333333333333333333333333333 
U 2 ,i = 0.10416666666666666666666666666666666667 
7 ^ 2,2 = —1.2291666666666666666666666666666666667 
TZ'2 } 3 = 1.1458333333333333333333333333333333333 
K 2 a = -0.020833333333333333333333333333333333333 

72 3 .1 = 0.020833333333333333333333333333333333333 
72 3 ’ 3 = -1.1458333333333333333333333333333333333 

72 3 . 4 = 1.2083333333333333333333333333333333333 
72 3 ’ 5 = -0.083333333333333333333333333333333333333 
72 4 ’ 3 = 0.083333333333333333333333333333333333333 
72 4 ’ 4 = -1.2500000000000000000000000000000000000 
72 4 ’ 5 = 1.2500000000000000000000000000000000000 
72 4 ’ 6 = -0.083333333333333333333333333333333333333 

The 7 Zi j * dx nonzero terms are: 

71.1.1 = -1.1250000000000000000000000000000000000 
72-1,2 = 1.2291666666666666666666666666666666667 
72. 4 ’ 3 = -0.083333333333333333333333333333333333333 

72 4 . 4 = -0.020833333333333333333333333333333333333 
72-2,1 = 1.2291666666666666666666666666666666667 
TZ-2,2 = -2.4583333333333333333333333333333333333 
72 - 2,3 = 1.2291666666666666666666666666666666667 
72 - 3,1 = -0.083333333333333333333333333333333333333 
72-3,2 = 1.2291666666666666666666666666666666667 
72 - 3,3 = —2.2916666666666666666666666666666666667 
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ft 3)4 = 1.2291666666666666666666666666666666667 
n 3 , 5 = -0.083333333333333333333333333333333333333 
K4,i = -0.020833333333333333333333333333333333333 
7Z 4) 3 = 1.2291666666666666666666666666666666667 
n 4A = -2.4583333333333333333333333333333333333 
TZ 4 , 5 = 1.3333333333333333333333333333333333333 
= -0.083333333333333333333333333333333333333 
n 5,3 = -0.083333333333333333333333333333333333333 
n 5 . 4 = 1.3333333333333333333333333333333333333 
TZ 5i5 = -2.5000000000000000000000000000000000000 
= 1.3333333333333333333333333333333333333 
= -0.083333333333333333333333333333333333333 


The BSij nonzero terms are: 

BSi,i = 1.8333333333333333333333333333333333333 
BSi, 2 = -3.0000000000000000000000000000000000000 
BS h 3 = 1.5000000000000000000000000000000000000 
BSi A = -0.33333333333333333333333333333333333333 


E.2.2 (3-6-3) Scheme Skew- Symmetric 

The flux point coefficients at point j are: 

Vj-2*Wj-2 = 0.01666666666666666666666666667*a+0. 0000000000000000000000000000 
Vj-2*Wj-i = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
Vj-2*Wj + o = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
Vj-2*Wj+i = — 0.01666666666666666666666666667*a+0. 01666666666666666666666666667 
Vj-2*Wj+2 = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
Vj-2*Wj+3 = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
2 = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
Vj—i*Wj— i = —0. 1333333333333333333333333333*a+0. 0000000000000000000000000000 
Vj-i*w j+ 0 = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
Vj-i*w j+4 = 0.1500000000000000000000000000*a-0.1500000000000000000000000000 
Vj-i*Wj+2 = — 0.01666666666666666666666666667*a+0. 01666666666666666666666666667 

Vj- i*w j+ 3 = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 

Vj+o*Wj-2 = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
V]+i)*w 3 -i = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
Vj + o*Wj + o = 0.6166666666666666666666666667*a+0. 0000000000000000000000000000 
v j+0 *w j+ i = —0.7500000000000000000000000000*a+0. 7500000000000000000000000000 
v j+ o*Wj+2 = 0.1500000000000000000000000000*a-0.1500000000000000000000000000 
Vj+Q*Wj+3 = — 0.01666666666666666666666666667*a+0. 01666666666666666666666666667 
Vj+i*Wj-2 = — 0.01666666666666666666666666667*a+0. 01666666666666666666666666667 
v j+ i *Wj-i = 0.1500000000000000000000000000*a-0.1500000000000000000000000000 
v j+1 *w j+0 = —0.7500000000000000000000000000*a+0. 7500000000000000000000000000 
Vj+i*Wj+i = 0.6166666666666666666666666667*a+0. 0000000000000000000000000000 
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v j+1 *w j+ 2 = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
Vj +1 *Wj+3 = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
Vj+2*Wj-2 = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
Vj+2*Wj-i = — 0.01666666666666666666666666667*a+0. 01666666666666666666666666667 
Vj+2*Wj +0 = 0.1500000000000000000000000000*a-0.1500000000000000000000000000 
v j+ 2 *w j+1 = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
Vj + 2*w j+ 2 = —0. 1333333333333333333333333333*a+0. 0000000000000000000000000000 
Vj + 2 *w j+ 3 = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
v j+3 *Wj-2 = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
Vj +i *Wj-i = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
Vj+$*Wj+Q = — 0.01666666666666666666666666667*a+0. 01666666666666666666666666667 
v j+3 *w j+1 = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
Vj+3*Wj+2 = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
Vj+3*Wj+s = 0.01666666666666666666666666667*a+0. 0000000000000000000000000000 


Flux point j=6 

v 3 *w 3 = —0.01396219135802469135802469136*a+0. 0000000000000000000000000000 
v 3 * w A = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 3 * W 5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 3 * w e = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

v 3 *w 7 = 0.01396219135802469135802469136*a— 0.01396219135802469135802469136 
v 3 *w 8 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 3 * wg = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

v 4 * w 3 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v A *w A = 0. 06997685185185185185185185185 *a+ 0.0000000000000000000000000000 
v A * W 5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v A * w 6 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 4 *w 7 = —0.06997685185185185185185185185*a+0. 06997685185185185185185185185 
v A * wg = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V A * wg = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 5 * xu 3 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 5 *w A = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 5 *u> 5 = —0.2056828703703703703703703704*a+0. 0000000000000000000000000000 

V 5 * w e = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

V 5 * w 7 = 0.2223495370370370370370370370 * a - 0.2223495370370370370370370370 
v 3 *w 3 = — 0.01666666666666666666666666667*a+0. 01666666666666666666666666667 
v 5 * W g = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 6 * xu 3 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 6 * w A = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 6 *UJ 5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 6 * W6 = 0.6496682098765432098765432099 * a + 0.0000000000000000000000000000 
v 6 *w 7 = —0.7830015432098765432098765432*a+0. 7830015432098765432098765432 
v 6 * W 8 = 0.1500000000000000000000000000 * a - 0.1500000000000000000000000000 
vgt- xjjg = — 0.01666666666666666666666666667*a+0. 01666666666666666666666666667 
v 7 *w 3 = 0.01396219135802469135802469136*a— 0.01396219135802469135802469136 
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v 7 *w A = —0.06997685185185185185185185185*a+0. 06997685185185185185185185185 
V7 * W5 = 0.2223495370370370370370370370 * a - 0.2223495370370370370370370370 
v 7 *w 6 = —0.7830015432098765432098765432*a+0. 7830015432098765432098765432 
v 7 *w 7 = 0.6166666666666666666666666667 * a + 0.0000000000000000000000000000 
v 7 *w 8 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 

V 7 * w 9 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

v s *w 3 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
v 8 * w A = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 8 *w§ = — 0.01666666666666666666666666667*a+0. 01666666666666666666666666667 
v 8 * w 6 = 0.1500000000000000000000000000 * a - 0.1500000000000000000000000000 
V S * w 7 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 8 *w 8 = —0. 1333333333333333333333333333 =t=a+0. 0000000000000000000000000000 

v 8 * wg = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
vg * '( 0 3 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
Vg * WA = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
Vg * W 5 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 

vg*we = — 0.01666666666666666666666666667*a+0. 01666666666666666666666666667 

Vg * w 7 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
Vg * w 8 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 

vg * wg = 0.01666666666666666666666666667 * a+ 0.0000000000000000000000000000 


Flux point j=5 

1,2*11, 2 = —0.02272762345679012345679012346*a+0. 0000000000000000000000000000 

1,2 * it , 3 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
1,2 * 11,4 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
1,2 * it, 5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

1,2*11,6 = 0.02272762345679012345679012346*a— 0.02272762345679012345679012346 

V 2 * W 7 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

V 2 *w 8 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v-i * 102 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
1,3*11,3 = 0.08048225308641975308641975309 *a+ 0.0000000000000000000000000000 
v :i * w 4 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 3 * 11 , 5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

1,3*11, 6 = —0.09444444444444444444444444444*a+0. 09444444444444444444444444444 
i,3*ii, 7 = 0.01396219135802469135802469136*a— 0.01396219135802469135802469136 
v 3 *w 8 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
174 * w 2 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
v A * 11,3 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

1,4*11,4 = —0. 1737924382716049382716049383 =t=a+0. 0000000000000000000000000000 
v A * w 5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
1,4 * 11, 6 = 0.2437692901234567901234567901 * a - 0.2437692901234567901234567901 
1,4*11,7 = —0.06997685185185185185185185185*a+0. 06997685185185185185185185185 

1,4 * u, 8 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

v 5 *w 2 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 

V 5 * 11,3 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
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v 5 * „, 4 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 5 * W5 = 0.6160378086419753086419753086 * a + 0.0000000000000000000000000000 
v 5 *w 6 = —0. 8217206790123456790123456790 =t=a+0. 8217206790123456790123456790 
v 5 * id j = 0.2223495370370370370370370370 * a - 0.2223495370370370370370370370 
v§*w 8 = — 0.01666666666666666666666666667*a+0. 01666666666666666666666666667 
v 6 *w 2 = 0.02272762345679012345679012346*a— 0.02272762345679012345679012346 
v e *w 3 = —0.09444444444444444444444444444*a+0. 09444444444444444444444444444 
v 6 * W4 = 0.2437692901234567901234567901 * a - 0.2437692901234567901234567901 
v 6 *w 5 = —0. 8217206790123456790123456790 =t=a+0. 8217206790123456790123456790 
v 6 * W6 = 0.6496682098765432098765432099 * a + 0.0000000000000000000000000000 
v 6 * W7 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 6 *w 8 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 7 * w 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 7 *w 3 = 0.01396219135802469135802469136*a— 0.01396219135802469135802469136 
v 7 *w 4 = —0.06997685185185185185185185185*a+0. 06997685185185185185185185185 
v 7 * W5 = 0.2223495370370370370370370370 * a - 0.2223495370370370370370370370 
V7 * W6 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 7 *w 7 = —0. 1663348765432098765432098765 =t=a+0. 0000000000000000000000000000 
v 7 * w 8 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

V 8 * W 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

v 8 * (03 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

V 8 * W 4 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

v 8 *w§ = — 0.01666666666666666666666666667*a+0. 01666666666666666666666666667 
v 8 * W6 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

V 8 * W7 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

v 8 *w 8 = 0.01666666666666666666666666667 * a+ 0.0000000000000000000000000000 


Flux point j=4 

vi*wi = 0.06527777777777777777777777778 * a+ 0.0000000000000000000000000000 
Vi *w 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V) *w 3 = 0.0000000000000000000000000000* a + 0.0000000000000000000000000000 
V! *w A = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 

v 7 *w 5 = —0.06527777777777777777777777778*a+0. 06527777777777777777777777778 
V) *w 6 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
V! *w 7 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
v 2 * Wi = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 2 *w 2 = —0.1170100308641975308641975309*a+0. 0000000000000000000000000000 
v 2 * w :i = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 2 * w 4 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 2 *w 5 = 0.09428240740740740740740740741*a-0. 09428240740740740740740740741 
v 2 *w 6 = 0.02272762345679012345679012346*a— 0.02272762345679012345679012346 
v 2 * iu 7 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 3 * wi = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 3 * w 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 3 *w 3 = —0.1723919753086419753086419753*a+0. 0000000000000000000000000000 
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v 3 * Wa = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 3 * W 5 = 0.2528742283950617283950617284 * a - 0.2528742283950617283950617284 
v 3 *w 6 = —0.09444444444444444444444444444*a+0. 09444444444444444444444444444 
v 3 *w 7 = 0.01396219135802469135802469136*a— 0.01396219135802469135802469136 
v A * Wx = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v A * w 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 4 * w 3 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v A * w A = 0.7241242283950617283950617284 * a + 0.0000000000000000000000000000 
v A *w 3 = — 0.8979166666666666666666666667*a+0. 8979166666666666666666666667 
VA * W6 = 0.2437692901234567901234567901 * a - 0.2437692901234567901234567901 
v A *w 7 = —0.06997685185185185185185185185*a+0. 06997685185185185185185185185 
v 5 *w A = —0.06527777777777777777777777778*a+0. 06527777777777777777777777778 
v 5 *w 2 = 0.09428240740740740740740740741*a-0. 09428240740740740740740740741 
v 5 * w 3 = 0.2528742283950617283950617284 * a - 0.2528742283950617283950617284 
vfj * w 4 = —0.8979166666666666666666666667*a+0. 8979166666666666666666666667 
v 5 * W5 = 0.6160378086419753086419753086 * a + 0.0000000000000000000000000000 
v 5 * We = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 5 * w 7 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 6 * Wx = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 6 *w 2 = 0.02272762345679012345679012346*a— 0.02272762345679012345679012346 
v 6 *w 3 = —0.09444444444444444444444444444*a+0. 09444444444444444444444444444 
v 6 *w A = 0.2437692901234567901234567901 * a - 0.2437692901234567901234567901 
v 6 * W5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 6 *w 6 = —0. 1720524691358024691358024691 *a+0. 0000000000000000000000000000 
v 6 *w 7 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 7 * Wx = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 7 * w 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 7 *w 3 = 0.01396219135802469135802469136*a— 0.01396219135802469135802469136 
v 7 *w A = —0.06997685185185185185185185185*a+0. 06997685185185185185185185185 
v 7 *w 5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 7 * Wfi = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 7 *w 7 = 0.05601466049382716049382716049 * a+ 0.0000000000000000000000000000 


Flux point j=3 

V!*wi = —0.09305555555555555555555555556*a+0. 0000000000000000000000000000 
V! *w 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V! *w 3 = 0.0000000000000000000000000000* a + 0.0000000000000000000000000000 
Vl * w A = 0.1583333333333333333333333333 * a - 0.1583333333333333333333333333 
v A *w 5 = —0.06527777777777777777777777778*a+0. 06527777777777777777777777778 
V! *w 6 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 

V! *w 7 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
v 2 *wi = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 2 * w 2 = 0.2480478395061728395061728395 * a + 0.0000000000000000000000000000 
V 2 * (03 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 2 *w A = —0.3650578703703703703703703704*a+0. 3650578703703703703703703704 
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v 2 *w 5 = 0.09428240740740740740740740741*a-0. 09428240740740740740740740741 
v 2 *w 6 = 0.02272762345679012345679012346*a— 0.02272762345679012345679012346 
v 2 * w 7 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 3 * wi = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 3 * w 2 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
v 3 * w 3 = 0.3450077160493827160493827161 * a + 0.0000000000000000000000000000 
v 3 *w A = —0.5173996913580246913580246914*a+0. 5173996913580246913580246914 
v 3 *w 5 = 0.2528742283950617283950617284 * a - 0.2528742283950617283950617284 
v 3 *w e = —0.09444444444444444444444444444*a+0. 09444444444444444444444444444 
v 3 *w 7 = 0.01396219135802469135802469136*a— 0.01396219135802469135802469136 
v A * Wl = 0.1583333333333333333333333333 * a - 0.1583333333333333333333333333 
v A *w 2 = —0.3650578703703703703703703704*a+0. 3650578703703703703703703704 
v A *w 3 = —0.5173996913580246913580246914*a+0. 5173996913580246913580246914 
v A * WA = 0.7241242283950617283950617284 * a + 0.0000000000000000000000000000 
v A *w 5 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 

V A * wq = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

V 4 *w 7 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

v 5 *w A = —0.06527777777777777777777777778*a+0. 06527777777777777777777777778 
v 5 *w 2 = 0.09428240740740740740740740741*a-0. 09428240740740740740740740741 
v 5 * W3 = 0.2528742283950617283950617284 * a - 0.2528742283950617283950617284 
v 5 * WA = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 5 *w 5 = —0. 2818788580246913580246913580 =t=a+0. 0000000000000000000000000000 

V 5 * w e = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 5 *w 7 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 6 *w 1 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
v 6 *w 2 = 0.02272762345679012345679012346*a— 0.02272762345679012345679012346 
v 6 *w 3 = —0.09444444444444444444444444444*a+0. 09444444444444444444444444444 
v e * WA = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
V 6 * w 5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

v 6 *w 6 = 0.07171682098765432098765432099 * a+ 0.0000000000000000000000000000 
v 6 * w 7 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 7 * Wl = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 7 * w 2 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
v 7 *w 3 = 0.01396219135802469135802469136*a— 0.01396219135802469135802469136 
v 7 * w A = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

V 7 * u, 5 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
V 7 * w 6 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

v 7 *w 7 = —0.01396219135802469135802469136*a+0. 0000000000000000000000000000 


Flux point j=2 

vi*wi = —0.1638888888888888888888888889*a+0. 0000000000000000000000000000 
v A *w 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
vi*w 3 = 0.07083333333333333333333333333*a— 0.07083333333333333333333333333 
v A * w A = 0.1583333333333333333333333333 * a - 0.1583333333333333333333333333 
va*w 5 = —0.06527777777777777777777777778*a+0. 06527777777777777777777777778 
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V! *w 6 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 2 * m = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 2 *w 2 = 0.6638888888888888888888888889 * a + 0.0000000000000000000000000000 
v 2 *w 3 = -0.4158410493827160493827160494*a+0.4158410493827160493827160494 
v 2 *w 4 = —0. 3650578703703703703703703704 =t=a+0. 3650578703703703703703703704 
v 2 *w 5 = 0.09428240740740740740740740741*a-0. 09428240740740740740740740741 
v 2 *w 6 = 0.02272762345679012345679012346*a— 0.02272762345679012345679012346 
v ^ Wl = 0. 07083333333333333333333333333*a— 0.07083333333333333333333333333 
v 3 *w 2 = — 0.4158410493827160493827160494*o+0.4158410493827160493827160494 
v 3 * - (03 = 0.3450077160493827160493827161 * o + 0.0000000000000000000000000000 
v 3 * W4 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 3 * u, 5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 3 * wq = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V A * Wl = 0.1583333333333333333333333333 * a - 0.1583333333333333333333333333 
v A *w 2 = —0.3650578703703703703703703704*a+0. 3650578703703703703703703704 
v 4 *w 3 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 4 *w 4 = 0.2067245370370370370370370370 * a + 0.0000000000000000000000000000 
V 4 * W 5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 4 * w e = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v b *wi = —0.06527777777777777777777777778*a+0. 06527777777777777777777777778 
v 5 *w 2 = 0.09428240740740740740740740741*a-0. 09428240740740740740740740741 
v 5 * w 3 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 5 * W 4 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 5 *w 5 = —0.02900462962962962962962962963*a+0. 0000000000000000000000000000 
v 5 * w 6 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 6 * Wl = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 6 *ui 2 = 0.02272762345679012345679012346*a— 0.02272762345679012345679012346 
v 6 * W:i = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 6 * w 4 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

V 6 * W 5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

V(i*vj(j = —0.02272762345679012345679012346*Q!+0. 0000000000000000000000000000 


Flux point j=l 

V!*W! = 0.5000000000000000000000000000* a + 0.0000000000000000000000000000 
Vi *102 = —0.6638888888888888888888888889*a+0. 6638888888888888888888888889 
Vl *w 3 = 0.07083333333333333333333333333*a— 0.07083333333333333333333333333 
Vi * w 4 = 0.1583333333333333333333333333 * a - 0.1583333333333333333333333333 
Vi*to 5 = —0.06527777777777777777777777778*a+0. 06527777777777777777777777778 
v 2 * Wl = —0.6638888888888888888888888889*a+0. 6638888888888888888888888889 
v 2 *w 2 = 0.6638888888888888888888888889 * a + 0.0000000000000000000000000000 

V 2 *w 3 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 2 *W 4 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 2 * w 5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

v 3 *«7i = 0.07083333333333333333333333333*a— 0.07083333333333333333333333333 
v 3 *W 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
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v 3 *w 3 = —0.07083333333333333333333333333*a+0. 0000000000000000000000000000 
v 3 * w A = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 3 * u, 5 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
v A * m = 0.1583333333333333333333333333 * a - 0.1583333333333333333333333333 
v A * w 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v A * w 3 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v A *w A = —0. 1583333333333333333333333333 =t=a+0. 0000000000000000000000000000 
v A *w 5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 5 * Wl = —0.06527777777777777777777777778*a+0. 06527777777777777777777777778 
v 5 * w 2 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
v 5 * w 3 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
v 5 * WA = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 5 *w 5 = 0.06527777777777777777777777778 * a+ 0.0000000000000000000000000000 


Flux point j=0 

Vi*wi = 0.0000000000000000000000000000 * a + 1.000000000000000000000000000 


The P-norm V is: 

Vij/dx = 0.30833333333333333333333333333333333333 
V 2 ,2 /dx = 1.4360879629629629629629629629629629630 
V 3 \ 3 /dx = 0.51331018518518518518518518518518518519 
V AA /dx = 1.3928240740740740740740740740740740741 
Vst/dx = 0.79745370370370370370370370370370370370 
Vf-, 6 /dx = 1.0596064814814814814814814814814814815 
'Prj/dx = 0.99238425925925925925925925925925925926 
V 8 ,8 /dx = 1.0000000000000000000000000000000000000 

The P-norm inverse V is: 

V x \ * dx = 3.2432432432432432432432432432432432432 

* dx = 0.69633617563145763148986927577813955737 

* dx = 1.9481397970687711386696730552423900789 
V A l * dx = 0.71796576366960279208908093734419145754 

* dx = 1.2539912917271407837445573294629898403 

* dx = 0.94374658656471873293282359366466411797 
Vjj * dx = 1.0076741853467378880828532107951762264 

* dx = 1.0000000000000000000000000000000000000 


The Qij * dx nonzero terms are: 

Qi .1 = -0.50000000000000000000000000000000000000 
Qi,2 = 0.66388888888888888888888888888888888889 
Qi’ 3 = -0.070833333333333333333333333333333333333 
Qi , 4 = -0.15833333333333333333333333333333333333 
Qi’ 5 = 0.065277777777777777777777777777777777778 
Q 2 i = -0.66388888888888888888888888888888888889 
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Q 2 , 3 = 0.41584104938271604938271604938271604938 
0,2,4: = 0.36505787037037037037037037037037037037 
Q 2 ’ 5 = -0.094282407407407407407407407407407407407 
Q 2 ’ 6 = -0.022727623456790123456790123456790123457 
Q 3 , i = 0.070833333333333333333333333333333333333 
Q 3 ,2 = -0.41584104938271604938271604938271604938 

0 3 .4 = 0.51739969135802469135802469135802469136 

0 3 . 5 = -0.25287422839506172839506172839506172840 
Q 3 ,6 = 0.094444444444444444444444444444444444444 
Q 3 j = -0.013962191358024691358024691358024691358 
Q 4 ’i = 0.15833333333333333333333333333333333333 
Q 4 ’ 2 = -0.36505787037037037037037037037037037037 
Q 4 ’ 3 = -0.51739969135802469135802469135802469136 
Q\, 5 = 0.89791666666666666666666666666666666667 

04.6 = -0.24376929012345679012345679012345679012 

04.7 = 0.069976851851851851851851851851851851852 

05,4 = -0.065277777777777777777777777777777777778 
Q 5 ’ 2 = 0.094282407407407407407407407407407407407 

05.3 = 0.25287422839506172839506172839506172840 

05.4 = —0.89791666666666666666666666666666666667 
Q 5 ,6 = 0.82172067901234567901234567901234567901 
Q 5 ,7 = -0.22234953703703703703703703703703703704 

05.8 = 0.016666666666666666666666666666666666667 

06.2 = 0.022727623456790123456790123456790123457 
Q 6 ,3 = -0.094444444444444444444444444444444444444 

06.4 = 0.24376929012345679012345679012345679012 
Q 6 , 5 = -0.82172067901234567901234567901234567901 
06,7 = 0.78300154320987654320987654320987654321 
g 6 ’ 8 = -0.15000000000000000000000000000000000000 

0 6 .9 = 0.016666666666666666666666666666666666667 

07.3 = 0.013962191358024691358024691358024691358 

07.4 = -0.069976851851851851851851851851851851852 

07. 5 = 0.22234953703703703703703703703703703704 

07.6 = -0.78300154320987654320987654320987654321 
Q 7 ’ 8 = 0.75000000000000000000000000000000000000 
Q 7 ’ 9 = -0.15000000000000000000000000000000000000 

07.10 = 0.016666666666666666666666666666666666667 
Qs, 5 = —0.016666666666666666666666666666666666667 
0 8 , e = 0.15000000000000000000000000000000000000 
Os,7 = -0.75000000000000000000000000000000000000 
Os,9 = 0.75000000000000000000000000000000000000 
Q 8 ’ 10 = -0.15000000000000000000000000000000000000 
Qs t ii = 0.016666666666666666666666666666666666667 

The (3-6-3) viscous terms. 
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The IZij * dx nonzero terms: 

7*0,1 = -2.2833333333333333333333333333333333333 
n 02 = 5.0000000000000000000000000000000000000 
7* 0 ’ 3 = -5.0000000000000000000000000000000000000 
Koa = 3.3333333333333333333333333333333333333 
7* 0 ’ 5 = -1.2500000000000000000000000000000000000 
Kofi = 0.20000000000000000000000000000000000000 
7* r’l = -1.1840277777777777777777777777777777778 
7* i’ 2 = 1.3277777777777777777777777777777777778 
7* i’ 3 = -0.070833333333333333333333333333333333333 
7* M = -0.10555555555555555555555555555555555556 
7* i’ 5 = 0.032638888888888888888888888888888888889 
TZ - 2,1 = 0.14375000000000000000000000000000000000 
n-2,2 = -1.1225212191358024691358024691358024691 
7* 2 ’ 3 = 0.76084876543209876543209876543209876543 
U 2 a = 0.25950231481481481481481481481481481481 
7* 2 ’ 5 = -0.030216049382716049382716049382716049383 
7* 2 ’ 6 = -0.011363811728395061728395061728395061728 
7* 3,1 = 0.072916666666666666666666666666666666667 
K : i . 2 = -0.29083912037037037037037037037037037037 
n 3 A = -0.83790702160493827160493827160493827160 
n 3A = 1.2943016975308641975308641975308641975 
TZ.iA = -0.28309027777777777777777777777777777778 
7* 3 ’ 6 = 0.051599151234567901234567901234567901235 
1 l 3 j = -0.0069810956790123456790123456790123456790 
7 * 4 ’i = -0.032638888888888888888888888888888888889 
7* 4 ’ 2 = 0.074218750000000000000000000000000000000 
7* 4 ’ 3 = 0.19689236111111111111111111111111111111 
K 4 a = -1.5987152777777777777777777777777777778 
n 4 A = 1.5127430555555555555555555555555555556 
7* 4j6 = -0.19217013888888888888888888888888888889 
7* 4 ’ 7 = 0.039670138888888888888888888888888888889 
7* 5 ’ 2 = 0.011363811728395061728395061728395061728 
7* 5 ’ 3 = -0.055981867283950617283950617283950617284 
7*5,4 = 0.19711805555555555555555555555555555556 
7* 5 ’ 5 = -1.4322029320987654320987654320987654321 
7* 5 ’ 6 = 1.4512712191358024691358024691358024691 
7*57 = -0.18267939814814814814814814814814814815 
7 * 5^8 = 0.011111111111111111111111111111111111111 
7* 6 ’ 3 = 0.0069810956790123456790123456790123456790 
7 * m = -0.046651234567901234567901234567901234568 
7* 6 ’ 5 = 0.21123842592592592592592592592592592593 
7* 6 ’ 6 = -1.4271141975308641975308641975308641975 
7 * 6 j = 1.3833236882716049382716049382716049383 
7* 6 ’ 8 = -0.13888888888888888888888888888888888889 
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ft 6 ,g = 0.011111111111111111111111111111111111111 

n 7A = -0.011111111111111111111111111111111111111 

n 7 ]e = 0.13888888888888888888888888888888888889 
1Z 7 j = -1.3611111111111111111111111111111111111 
H 7 [s = 1.3611111111111111111111111111111111111 
n 7 ,g = -0.13888888888888888888888888888888888889 
H 7 ]io = 0.011111111111111111111111111111111111111 


The IZij * dx nonzero terms: 

K lA = -1.1840277777777777777777777777777777778 
Ki i2 = 1.3277777777777777777777777777777777778 
nl 3 = -0.070833333333333333333333333333333333333 
Ui A = -0.10555555555555555555555555555555555556 
Hi A = 0.032638888888888888888888888888888888889 
H 2 ,i = 1.3277777777777777777777777777777777778 
U 2 ,2 = -2.4502989969135802469135802469135802469 
TZ 2 ,3 = 0.83168209876543209876543209876543209877 
U 2A = 0.36505787037037037037037037037037037037 
7Z 2}5 = -0.062854938271604938271604938271604938272 
n 2fi = -0.011363811728395061728395061728395061728 
77,34 = -0.070833333333333333333333333333333333333 
^ 3 ’ 2 = 0.83168209876543209876543209876543209877 
n 3A = -1.5987557870370370370370370370370370370 
77.,3, 4 = 1.0347993827160493827160493827160493827 
77 3 ’ 5 = -0.25287422839506172839506172839506172840 
77 3 ’ 6 = 0.062962962962962962962962962962962962963 
T7 3 ’ 7 = -0.0069810956790123456790123456790123456790 
77 41 = -0.10555555555555555555555555555555555556 
77 4 ’ 2 = 0.36505787037037037037037037037037037037 
77 4 ’ 3 = 1.0347993827160493827160493827160493827 
n 4A = -2.8930169753086419753086419753086419753 
TZ 4 ,5 = 1.7958333333333333333333333333333333333 
= -0.24376929012345679012345679012345679012 
= 0.046651234567901234567901234567901234568 
U 5A = 0.032638888888888888888888888888888888889 
n 5 . 2 = -0.062854938271604938271604938271604938272 
U 5A = -0.25287422839506172839506172839506172840 
n 5A = 1.7958333333333333333333333333333333333 
n 5A = -2.9449459876543209876543209876543209877 
= 1.6434413580246913580246913580246913580 
TZ 5 j = -0.22234953703703703703703703703703703704 

n 5 . 8 = 0.011111111111111111111111111111111111111 

TZ 6:2 = -0.011363811728395061728395061728395061728 
n 6:3 = 0.062962962962962962962962962962962962963 
U 6A = -0.24376929012345679012345679012345679012 
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^ 6j5 = 1.6434413580246913580246913580246913580 
= -2.8783854166666666666666666666666666667 
n 6:7 = 1.5660030864197530864197530864197530864 
= -0.15000000000000000000000000000000000000 

n^ 9 = o.oiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiin 

U 7 [ 3 = -0.0069810956790123456790123456790123456790 
TZ 7A = 0.046651234567901234567901234567901234568 
n 7 [ 5 = -0.22234953703703703703703703703703703704 
n 7fi = 1.5660030864197530864197530864197530864 
n 7 [ 7 = -2.7444347993827160493827160493827160494 
n 7 [ 8 = 1.5000000000000000000000000000000000000 
n 7 [g = -0.15000000000000000000000000000000000000 
Hr ,10 = 0.011111111111111111111111111111111111111 

n 8 ,5 = o.oiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiin 

n 8 ,6 = -0.15000000000000000000000000000000000000 
n 8:7 = 1.5000000000000000000000000000000000000 
n 8t8 = -2.7222222222222222222222222222222222222 
^8,9 = 1.5000000000000000000000000000000000000 
ft 8 ’io = -0.15000000000000000000000000000000000000 

KS'U = 0.011111111111111111111111111111111111111 


The BSi'j nonzero terms 

BSi'i = 2.2833333333333333333333333333333333333 
BS h 2 = -5.0000000000000000000000000000000000000 
BSi ,3 = 5.0000000000000000000000000000000000000 
BSi A = -3.3333333333333333333333333333333333333 
BSi'O = 1.2500000000000000000000000000000000000 
BSi'O = -0.20000000000000000000000000000000000000 


E.2.3 (4-8-4) Scheme Skew- Symmetric 

The coefficients for flux point j are : 

Vj-o*Wj-z = —0.003571428571428571428571428571*a+0. 0000000000000000000000000000 
Vj-^Wj-2 = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
Vj-3*Wj-i = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
Vj-3*Wj+o = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
Vj-3*w j+ i = 0.003571428571428571428571428571*a— 0.003571428571428571428571428571 
Vj-3*Wj+2 = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
Vj-3*Wj+3 = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 

Vj-3* Wj +A = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 

Vj-2*Wj- 3 = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
Vj-2*Wj-2 = 0.03452380952380952380952380952*a:+0. 0000000000000000000000000000 
Vj-2*Wj-i = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
Vj-2*Wj+o = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
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Vj- 2 *Wj + i 
Vj-2*Wj+2 
Vj-2*Wj +3 
Vj-2*Wj+A 
Vj- l*Wj-3 
Vj-i*Wj- 2 
Vj-l*Wj-l 
Vj-i*W j+ o 
Vj-l*Wj + l 
Vj-l*W j+2 
Vj-i*W j+3 
Vj-i*W j+ 4 
Vj +0 *Wj-3 
V j+0 *Wj-2 
Vj+o*Wj-l 

Vj +0 *Wj + 0 

Vj+0*Wj+l 

Vj+0*Wj+2 

Vj+0*Wj + 3 
v j+0* w j+A 
V j+ i*Wj-3 
Vj + i*Wj- 2 
Vj+l*Wj-l 
Vj+l*Wj +0 
Vj + l*W j + l 
Vj+l*Wj+2 
Vj+l*W j+ 3 
Vj+l* w j+A 

V j+2 *Wj-3 

V j+2 *Wj-2 

Vj + 2*Wj-l 

Vj + 2*Wj + o 

Vj + 2*Wj + i 

Vj + 2*Wj + 2 

Vj + 2*Wj + 3 

v j+2* w j+4 

Vj+3*Wj-3 

V j+ 3*Wj-2 

V j+ 3*Wj~i 

v j+3* w j+0 

Vj + 3*Wj +1 

Vj+3*W j+ 2 

Vj + 3*Wj +3 

Vj+3*Wj + 4 

Vj+4*1Vj—3 


— 0. 0380952380952380952380952381 0*ct+0. 038095238095238095238095238 10 

0.003571428571428571428571428571*a— 0.003571428571428571428571428571 

0.0000000000000000000000000000*a+0.0000000000000000000000000000 

0.0000000000000000000000000000*a+0.0000000000000000000000000000 

0.0000000000000000000000000000*a+0.0000000000000000000000000000 

0.0000000000000000000000000000*a+0.0000000000000000000000000000 

— 0.1654761904761904761904761905*cH-0.0000000000000000000000000000 

0.0000000000000000000000000000*a+0.0000000000000000000000000000 

0.2000000000000000000000000000*a-0.2000000000000000000000000000 

—0.03809523809523809523809523810*a+0. 03809523809523809523809523810 

0.003571428571428571428571428571*a— 0.003571428571428571428571428571 

0.0000000000000000000000000000*a+0.0000000000000000000000000000 

0.0000000000000000000000000000*a+0.0000000000000000000000000000 

0.0000000000000000000000000000*a+0.0000000000000000000000000000 

0.0000000000000000000000000000*a+0.0000000000000000000000000000 

0.6345238095238095238095238095*a+0. 0000000000000000000000000000 

-0.8000000000000000000000000000*a+0.8000000000000000000000000000 

0.2000000000000000000000000000*a-0.2000000000000000000000000000 

—0.03809523809523809523809523810*a+0. 03809523809523809523809523810 

0.003571428571428571428571428571*a— 0.003571428571428571428571428571 

0.003571428571428571428571428571*a— 0.003571428571428571428571428571 

—0.03809523809523809523809523810*a+0. 03809523809523809523809523810 

0.2000000000000000000000000000*a-0.2000000000000000000000000000 

-0.8000000000000000000000000000*a+0.8000000000000000000000000000 

0.6345238095238095238095238095*a+0. 0000000000000000000000000000 

0.0000000000000000000000000000*a+0.0000000000000000000000000000 

0.0000000000000000000000000000*a+0.0000000000000000000000000000 

0.0000000000000000000000000000*a+0.0000000000000000000000000000 

0.0000000000000000000000000000*a+0.0000000000000000000000000000 

0.003571428571428571428571428571*a— 0.003571428571428571428571428571 

—0.03809523809523809523809523810*a+0. 03809523809523809523809523810 

0.2000000000000000000000000000*a-0.2000000000000000000000000000 

0.0000000000000000000000000000*a+0.0000000000000000000000000000 

— 0.1654761904761904761904761905*a+0.0000000000000000000000000000 

0.0000000000000000000000000000*a+0.0000000000000000000000000000 

0.0000000000000000000000000000*a+0.0000000000000000000000000000 

0.0000000000000000000000000000*a+0.0000000000000000000000000000 

0.0000000000000000000000000000*a+0.0000000000000000000000000000 

0.003571428571428571428571428571*a-0. 003571428571428571428571428571 

—0.03809523809523809523809523810*a+0. 03809523809523809523809523810 

0.0000000000000000000000000000*a+0.0000000000000000000000000000 

0.0000000000000000000000000000*a+0.0000000000000000000000000000 

0.03452380952380952380952380952*a+0. 0000000000000000000000000000 

0.0000000000000000000000000000*a+0.0000000000000000000000000000 

0.0000000000000000000000000000*a+0.0000000000000000000000000000 
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v j+ 4 *Wj -2 = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
Vj+4*Wj-i = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
v j+ 4 *w j+0 = 0.003571428571428571428571428571*a— 0.003571428571428571428571428571 
v j+4 *w j+ i = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
Vj + 4 *w j+2 = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
v j+ 4 *w j+3 = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
v j+ 4 *Wj +4 = —0.003571428571428571428571428571*a+0. 0000000000000000000000000000 


Flux point j=8 

v A *w A = —0. 0096027579627530024355421 18082*cH-0. 0000000000000000000000000000 

V A * W 5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V A * w e = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

v A *w 7 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v A *w 8 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v A *w 9 = 0. 009602757962753002435542118082*0;— 0.009602757962753002435542118082 

V A * Wl0 = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 

v A *w n = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 

V A * Wl2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

v 5 * w A = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 5 *w 5 = 0.04329003028680608045687410767* a+ 0.0000000000000000000000000000 
v 5 * m = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 5 * w 7 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 5 * w 8 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 5 *w 9 = —0.04329003028680608045687410767*a+0. 04329003028680608045687410767 
v 5 * Wl0 = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
v 5 *wu = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
v 5 * Wl2 = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
V 6 * w 4 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 6 * w 5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

v 6 *w 6 = —0.05535800500755857898715041572*a+0. 0000000000000000000000000000 
v 6 * w 7 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 6 *w 8 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 6 *u> 9 = 0.05178657643613000755857898715*a— 0.05178657643613000755857898715 
v 6 *w 10 = 0.003571428571428571428571428571*a— 0.003571428571428571428571428571 
v 6 * wu = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 

V 6 * Wl2 = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
V 7 * WA = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

V 7 * w 5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 7 * W6 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 7 *w 7 = —0.08327658231502477534223565970*a+0. 0000000000000000000000000000 
v 7 * w 8 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 7 * w 9 = 0.1178003918388342991517594692 * a - 0.1178003918388342991517594692 
v 7 *w 10 = —0.03809523809523809523809523810*a+0. 03809523809523809523809523810 
v 7 *w u = 0.003571428571428571428571428571*a— 0.003571428571428571428571428571 

V 7 * Wl2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
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v 8 * W 4 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 8 * W5 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
V 8 * vj (i = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 

V 8 * id j = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

v 8 *w 8 = 0.6049473149985302763080540858 * a + 0.0000000000000000000000000000 
v 8 *w 9 = —0. 7704235054747207524985302763 =t=a+0. 7704235054747207524985302763 
v 8 * w 10 = 0.2000000000000000000000000000 * a - 0.2000000000000000000000000000 
v 8 *wn = — 0.03809523809523809523809523810*a+0.03809523809523809523809523810 
v 8 *w 12 = 0.003571428571428571428571428571*a— 0.003571428571428571428571428571 
v 9 *w A = 0.009602757962753002435542118082*a— 0.009602757962753002435542118082 
v 9 *w 5 = -0.04329003028680608045687410767*a+0.04329003028680608045687410767 
v 9 *w 6 = 0.05178657643613000755857898715*a— 0.05178657643613000755857898715 
vg * w 7 = 0.1178003918388342991517594692 * a - 0.1178003918388342991517594692 
vg*w 8 = —0. 7704235054747207524985302763 =t=a+0. 7704235054747207524985302763 
V9 * W9 = 0.6345238095238095238095238095 * a + 0.0000000000000000000000000000 

v 9 * w 10 = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
Vg * w u = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 

Vg * W12 = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 

u 10 * w 4 = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
u 10 * io 5 = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
v 10 *iu 6 = 0.003571428571428571428571428571*0!— 0.003571428571428571428571428571 
= — 0.03809523809523809523809523810*a+0.03809523809523809523809523810 
u 10 *w 8 = 0.2000000000000000000000000000 * a - 0.2000000000000000000000000000 
vio *wg = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
v 10 *w 10 = —0.1654761904761904761904761905*a+0. 0000000000000000000000000000 
v 10 *wu = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
v w *w 12 = 0.0000000000000000000000000000*a+0. 0000000000000000000000000000 
Vll *iu 4 = 0.0000000000000000000000000000*a+ 0.0000000000000000000000000000 
vn*w 5 = 0.0000000000000000000000000000*a+ 0.0000000000000000000000000000 
v u *w 6 = 0.0000000000000000000000000000*a+ 0.0000000000000000000000000000 
vn *w 7 = 0.003571428571428571428571428571*a— 0.003571428571428571428571428571 
Vll *w 8 = —0.03809523809523809523809523810*a+0. 03809523809523809523809523810 
v n *w 9 = 0.0000000000000000000000000000*a+ 0.0000000000000000000000000000 
u n *«; 10 = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
vn*w u = 0.03452380952380952380952380952*a+0. 0000000000000000000000000000 
v u *wi2 = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
V12 * 104 = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
V 12 * w 5 = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
V 12 *W 6 = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 

V12 *W 7 = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 

v 12 *w 8 = 0.003571428571428571428571428571*a— 0.003571428571428571428571428571 
V12 *wg = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
Vi 2 *w w = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
V 12 *w u = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
v 12 *w 12 = —0.003571428571428571428571428571*a+0. 0000000000000000000000000000 
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Flux point j=7 

v 3 *w 3 = —0.09050595238095238095238095238*a+0. 0000000000000000000000000000 
v 3 * W4 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
v 3 * w 5 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 

V 3 * u, 6 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

^3 * w 7 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 3 *w 8 = 0.09050595238095238095238095238*a— 0.09050595238095238095238095238 

V 3 * „, 9 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

v 3 * wio = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
v 3 * w u = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
v A * (03 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 4 * W 4 = 0.4166642632747543461829176115 * a + 0.0000000000000000000000000000 
„ 4 * w 5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v A * w 6 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
v A * w 7 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v A *w 8 = —0. 4262670212375073486184597296 =t=a+0.4262670212375073486184597296 
v A *w 9 = 0.009602757962753002435542118082*a— 0.009602757962753002435542118082 
v A * wio = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
v A *wn= 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
v 5 *w 3 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

V 5 * w A = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 

v 5 *w 5 = -0.7125127557581674645166708659*a+0.0000000000000000000000000000 

V 5 * w 6 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 5 * W7 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 

v 5 * w 8 = 0.7558027860449735449735449735 * a - 0.7558027860449735449735449735 
v 5 *w 9 = -0.04329003028680608045687410767*a+0.04329003028680608045687410767 
V 5 * w 10 = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
B5*(UU = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
ve*w 3 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 6 * w A = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 6 *w 5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 6 * w 6 = 0.4365587627026119089611153103 * a + 0.0000000000000000000000000000 

V 6 *w 7 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

v 6 *w 8 = —0. 4919167677101704879482657260 =t=a+0.4919167677101704879482657260 
vq*u)q = 0.05178657643613000755857898715*a— 0.05178657643613000755857898715 
v 6 *w 10 = 0.003571428571428571428571428571*a— 0.003571428571428571428571428571 

BB*(UU = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
V 7 * ,o 3 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 7 * w A = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 7 *w 5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 7 * wq = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
v 7 * w 7 = 0.4497956821617535903250188964 * a + 0.0000000000000000000000000000 
v 7 *w 8 = —0. 5330722644767783656672545561 *a+0. 5330722644767783656672545561 
V7 * W9 = 0.1178003918388342991517594692 * a - 0.1178003918388342991517594692 
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V 7 *w 10 = — 0.03809523809523809523809523810*cH-0.03809523809523809523809523810 
v 7 *w u = 0.003571428571428571428571428571*0!— 0.003571428571428571428571428571 
v 8 *w 3 = 0. 09050595238095238095238095238*O!— 0.09050595238095238095238095238 
v 8 *w A = — 0.4262670212375073486184597296*O!+0.4262670212375073486184597296 
v 8 * ti7 5 = 0.7558027860449735449735449735 * a - 0.7558027860449735449735449735 
v 8 *w 6 = — 0.4919167677101704879482657260*O!+0.4919167677101704879482657260 
v 8 *w 7 = —0. 5330722644767783656672545561 *Q!+0. 5330722644767783656672545561 
v 8 * w 8 = 0.6049473149985302763080540858 * a + 0.0000000000000000000000000000 
v 8 * w 9 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 8 * « 7 10 = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
v 8 * wu = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
v 9 * w :i = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 9 *w A = 0. 009602757962753002435542118082*O!— 0.009602757962753002435542118082 
v 9 *w 5 = —0.04329003028680608045687410767*o+0. 04329003028680608045687410767 
v 9 *w 6 = 0.05178657643613000755857898715*o— 0.05178657643613000755857898715 
v 9 * w 7 = 0.1178003918388342991517594692 * o - 0.1178003918388342991517594692 
vg * w 8 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
vg *w 9 = — 0.13589969595091 12286890064668*a+0. 0000000000000000000000000000 
Vg * Wl0 = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
Vg * wu = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
u 10 * w 3 = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
u 10 * iv A = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
u 10 * w 5 = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 

= 0.003571428571428571428571428571*0!— 0.003571428571428571428571428571 
v w *iu 7 = —0.03809523809523809523809523810*Q!+0. 03809523809523809523809523810 
uio * w 8 = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 

u 10 * Wg = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 

V W *IU 10 = 0.03452380952380952380952380952*O!+0. 0000000000000000000000000000 

vio*w n = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
v u *W 3 = 0.0000000000000000000000000000*a+ 0.0000000000000000000000000000 

vu *io 4 = 0.0000000000000000000000000000*a+ 0.0000000000000000000000000000 
vu *w 5 = 0.0000000000000000000000000000*a+ 0.0000000000000000000000000000 
vu *w 6 = 0.0000000000000000000000000000*a+ 0.0000000000000000000000000000 
vu*w 7 = 0.003571428571428571428571428571*0!— 0.003571428571428571428571428571 
VU *w 8 = 0.0000000000000000000000000000*Q!+ 0.0000000000000000000000000000 
vu*wg = 0.0000000000000000000000000000*a+ 0.0000000000000000000000000000 
vu*wio = 0.0000000000000000000000000000*a+0.0000000000000000000000000000 
vu*wu = —0.003571428571428571428571428571*a+0. 0000000000000000000000000000 


Flux point j=6 

v 2 *W 2 = —0. 21072090275153271 18501721676*O!+0. 0000000000000000000000000000 
v 2 * w :i = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

V 2 * W 4 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 2 * W 5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 2 * W6 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
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v 2 * w 7 = 0.2107209027515327118501721676 * a - 0.2107209027515327118501721676 
v 2 * w% = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

V 2 * w 9 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

«2 * Wl0 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 3 * w 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 3 *w 3 = 1.090508276381540270429159318 * a + 0.0000000000000000000000000000 
v 3 * w 4 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 3 *u> 5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 3 * we = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 3 * w 7 = -1.181014228762492651381540270 * a + 1.181014228762492651381540270 
v 3 *w 8 = 0.09050595238095238095238095238*a— 0.09050595238095238095238095238 
v 3 * Wg = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

V 3 * «, 10 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

v A *w 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v A * w 3 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v A * w 4 = -2.220792581963340891912320484 * a + 0.0000000000000000000000000000 
V A * w 5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V A * w 6 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v A *w 7 = 2.637456845238095238095238095 * a - 2.637456845238095238095238095 
v A *w 8 = —0. 4262670212375073486184597296 =t=a+0.4262670212375073486184597296 
v A *w 9 = 0.009602757962753002435542118082*0;— 0.009602757962753002435542118082 

v A * Wl0 = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
v 5 * w 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 5 * w 3 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 5 *u> A = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 5 *w 5 = 2.067331343301209372637944067 * a + 0.0000000000000000000000000000 
v 5 * w 6 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 5 *w 7 = -2.779844099059376837154614932 * a + 2.779844099059376837154614932 
v 5 * W8 = 0.7558027860449735449735449735 * a - 0.7558027860449735449735449735 
v b *w 9 = —0.04329003028680608045687410767*a+0. 04329003028680608045687410767 
v 5 * «, 10 = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
ve*w 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
ve * w 3 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
ve * w A = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 6 * W5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
ve*we = —0.2263261349678760393046107332*a+0. 0000000000000000000000000000 
v 6 * W7 = 0.6628848976704879482657260435 * a - 0.6628848976704879482657260435 
v 6 *w 8 = — 0.4919167677101704879482657260*a+0.4919167677101704879482657260 
v 6 *u> 9 = 0.05178657643613000755857898715*a— 0.05178657643613000755857898715 
D6 *^io = 0.003571428571428571428571428571*a— 0.003571428571428571428571428571 
v 7 * w 2 = 0.2107209027515327118501721676 * a - 0.2107209027515327118501721676 
v 7 *w 3 = -1.181014228762492651381540270 * a + 1.181014228762492651381540270 
v 7 *w A = 2.637456845238095238095238095 * a - 2.637456845238095238095238095 
v 7 * w 5 = -2.779844099059376837154614932 * a + 2.779844099059376837154614932 
V7 * W6 = 0.6628848976704879482657260435 * a - 0.6628848976704879482657260435 
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v 7 * w 7 = 0.4497956821617535903250188964 *a + 0.0000000000000000000000000000 
v 7 * w 8 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 

V 7 * w 9 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 7 * «, 10 = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 

v 8 * w 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 8 *w 3 = 0.09050595238095238095238095238*a— 0.09050595238095238095238095238 
v 8 *w A = — 0.4262670212375073486184597296*a+0.4262670212375073486184597296 
v 8 *w 5 = 0.7558027860449735449735449735 * a - 0.7558027860449735449735449735 
v 8 *w 6 = — 0.4919167677101704879482657260*a+0.4919167677101704879482657260 

V 8 * w 7 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

v 8 *w 8 = 0. 07187505052175191064079952969 *a+ 0.0000000000000000000000000000 
v 8 * W g = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

V 8 * «, 10 = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 

v 9 * w 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
Vg * w :i = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
Vg*w 4 = 0.009602757962753002435542118082*a— 0.009602757962753002435542118082 
V g*w 5 = —0.04329003028680608045687410767*a+0. 04329003028680608045687410767 
v 9 *w 6 = 0.05178657643613000755857898715*a— 0.05178657643613000755857898715 
Vg * W7 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

Vg * w 8 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

v 9 *w 9 = — 0.0180993041 1207692953724699756*a+0. 0000000000000000000000000000 

Vg * w 10 = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 

u 10 *w 2 = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
vio *w 3 = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
u 10 * W 4 = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
uio *w 5 = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
v 10 *w 6 = 0.003571428571428571428571428571*a— 0.003571428571428571428571428571 
u 10 *w 7 = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
u 10 *w 8 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
u 10 * lug = 0.0000000000000000000000000000 * a+ 0.0000000000000000000000000000 
'U 10 *i^ 10 = —0.003571428571428571428571428571*a+0. 0000000000000000000000000000 


Flux point j=5 

V!*W! = —0. 1239583333333333333333333333 =t=a+0. 0000000000000000000000000000 
vg *w 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
Vg *w 3 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v { *w A = 0.0000000000000000000000000000* a + 0.0000000000000000000000000000 
Vg *w 5 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 

Vg * we = 0.1239583333333333333333333333 * a - 0.1239583333333333333333333333 

Vg *w 7 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
Vg *w 8 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
Vg *wg = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
v 2 * m = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 2 *w 2 = 1.048880436009280255312001344 * a + 0.0000000000000000000000000000 
v 2 * w :i = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
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V2 * W 4 
V 2 * Wq 
V 2 * Wq 
V2 * W7 
V2*VJg 
V 2 * Wq 
V 3 * W\ 
V 3 *W 2 
V 3 *W 3 
V 3 * W4 
v 3 * wq 
v 3 * w 6 

V 3 * Wj 
V 3 *W 8 -- 
V 3 * WQ 
V4 * W\ 
V4 * W2 
V4 * W 3 
V4 * W4 
V4 * Wq 
V4 * Wq 
V4 * W7 
V 4 *W 8 : 
V4*Wg - 
Vq * W\ 
Vq * W 2 
Vq * W 3 
Vq * W4 
Vq*Wq 
Vq * Wq 
Vq * W7 
Vq * Wg 
Vq*Wq - 
Vq * W\ 
Vq * W 2 
Vq * W 3 
Vq > 1 = W4 
Vq * Wq 
Vq*Wq -- 
Vq * W7 
VQ*Wg 
Vq * Wq 
V7 * VJ\ 
V7 * W2 
V7 * W 3 


■ 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
: -1.259601338760812967162173511 * a + 1.259601338760812967162173511 
= 0.2107209027515327118501721676 * a - 0.2107209027515327118501721676 
= 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 >1= a + 0.0000000000000000000000000000 
-3.083680753049676660787771899 * a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 > 1 = a + 0.0000000000000000000000000000 
= 4.174189029431216931216931217 * a - 4.174189029431216931216931217 
: -1.181014228762492651381540270 * a + 1.181014228762492651381540270 
0.09050595238095238095238095238=m— 0.09050595238095238095238095238 
: 0.0000000000000000000000000000 > 1 = a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 >1= a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 >1= a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
= 3.940000223739187032837826489 * a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 > 1 = a + 0.0000000000000000000000000000 
: -6.160792805702527924750146972 * a + 6.160792805702527924750146972 
= 2.637456845238095238095238095 * a - 2.637456845238095238095238095 
— 0.4262670212375073486184597296*a+0.4262670212375073486184597296 
0. 009602757962753002435542118082*0!— 0.009602757962753002435542118082 
: 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
-1 .281241573365457294028722600 * a + 0.0000000000000000000000000000 
= 3.348572916666666666666666667 * a - 3.348572916666666666666666667 
: -2.779844099059376837154614932 * a + 2.779844099059376837154614932 
= 0.7558027860449735449735449735 * a - 0.7558027860449735449735449735 
—0.04329003028680608045687410767*a+0. 04329003028680608045687410767 
= 0.1239583333333333333333333333 * a - 0.1239583333333333333333333333 
= -1.259601338760812967162173511 * a + 1.259601338760812967162173511 
= 4.174189029431216931216931217 * a - 4.174189029431216931216931217 
: -6.160792805702527924750146972 * a + 6.160792805702527924750146972 
= 3.348572916666666666666666667 * a - 3.348572916666666666666666667 
—0.2263261349678760393046107332*a+0. 0000000000000000000000000000 
: 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
: 0.2107209027515327118501721676 * a - 0.2107209027515327118501721676 
= -1.181014228762492651381540270 * a + 1.181014228762492651381540270 
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v 7 * w 4 = 2.637456845238095238095238095 * a - 2.637456845238095238095238095 
v 7 * w 5 = - 2.779844099059376837154614932 * a + 2.779844099059376837154614932 
v 7 * W6 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 7 * w 7 = 1.112680579832241538590744940 *a + 0.0000000000000000000000000000 
v 7 * w 8 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 

V 7 * w 9 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

v 8 *wi = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 8 *w 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 8 *w 3 = 0.09050595238095238095238095238*a— 0.09050595238095238095238095238 
v 8 *w 4 = — 0.4262670212375073486184597296*a+0.4262670212375073486184597296 
v 8 * w 5 = 0.7558027860449735449735449735 * a - 0.7558027860449735449735449735 
v 8 * W6 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 8 * w 7 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 8 *w 8 = —0.4200417171884185773074661964*a+0. 0000000000000000000000000000 
V 8 * w 9 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 9 *wi = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 9 * w 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 9 * w 3 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 9 *w 4 = 0.009602757962753002435542118082*a— 0.009602757962753002435542118082 
v 9 *w 5 = -0.04329003028680608045687410767*a+0.04329003028680608045687410767 
Vg * We = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 9 * w 7 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 9 *w 8 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 9 * w 9 = 0.03368727232405307802133198959 * a+ 0.0000000000000000000000000000 

Flux point j=4 

Vi*wi = 0.4209583333333333333333333333* a + 0.0000000000000000000000000000 
v ! * ii)2 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
V! *w 3 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
Vi *W4 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 

vi*w§ = — 0.5449166666666666666666666667*a+0. 5449166666666666666666666667 
vi * w 6 = 0.1239583333333333333333333333 * a - 0.1239583333333333333333333333 
V! *w 7 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
V! *w 8 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
V! *w 9 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
v 2 * Wx = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 2 *w 2 = -1 .964242213745065927605610145 * a + 0.0000000000000000000000000000 

V 2 * w :i = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 2 * w 4 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

v 2 * W 5 = 3.013122649754346182917611489 * a - 3.013122649754346182917611489 
v 2 *Wd = -1.259601338760812967162173511 * a + 1.259601338760812967162173511 
V2 * W7 = 0.2107209027515327118501721676 * a - 0.2107209027515327118501721676 

V 2 *w 8 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 2 * w 9 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 3 * Wl = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
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V 3 *W 2 
V 3 * W 3 
V 3 * W4 
V 3 * WQ 
V 3 * Wq 
V 3 * W7 
V 3 *W 8 ~- 
V 3 * Wg 
V4 * W 1 
V4 * W2 
V4 * W 3 
V4 * IV 4 
V4 * WQ 
V4 * Wq 
V 4 * W7 
V4 *W 8 
V4*Wg - 
Vq*W\ - 
Vq * W2 
Vq = 1 = W 3 
Vq * W4 
Vq*Wq 
vq * wq 
vq * W7 
vq * w 8 
Vq * Wg 
Vq * W 1 
Vs * W2 
vq * w 3 

Vq * W4 
Vq * Wq 
Vq * Wq 
Vq * W7 
VQ*W 8 
Vq * Wg 
V7 * W 1 
V7 * W2 
V7 * W 3 
V7 * W4 
V7 * Wq 
V7 * Wq 
V7 *W7 
V7 * W 8 
V7 * Wg 
V 8 * W 1 


: 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
= 3.207549800668724279835390946 * a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 >1= a + 0.0000000000000000000000000000 
: -6.291230553718400940623162845 * a + 6.291230553718400940623162845 
= 4.174189029431216931216931217 * a - 4.174189029431216931216931217 
: -1.181014228762492651381540270 * a + 1.181014228762492651381540270 
0.09050595238095238095238095238=m— 0.09050595238095238095238095238 
: 0.0000000000000000000000000000 >1= a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
-1.164265920256991685563114135 * a + 0.0000000000000000000000000000 
= 5.104266143996178718400940623 * a - 5.104266143996178718400940623 
: -6.160792805702527924750146972 * a + 6.160792805702527924750146972 
= 2.637456845238095238095238095 * a - 2.637456845238095238095238095 
— 0.4262670212375073486184597296*a+0.4262670212375073486184597296 
0.009602757962753002435542118082*a— 0.009602757962753002435542118082 
— 0.5449166666666666666666666667*a+0. 5449166666666666666666666667 
= 3.013122649754346182917611489 * a - 3.013122649754346182917611489 
: -6.291230553718400940623162845 * a + 6.291230553718400940623162845 
= 5.104266143996178718400940623 * a - 5.104266143996178718400940623 
-1 .281241573365457294028722600 * a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 >1= a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 >1= a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
0.1239583333333333333333333333 * a - 0.1239583333333333333333333333 
= -1.259601338760812967162173511 * a + 1.259601338760812967162173511 
= 4.174189029431216931216931217 * a - 4.174189029431216931216931217 
: -6.160792805702527924750146972 * a + 6.160792805702527924750146972 
: 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
= 3.122246781698790627362055934 * a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 > 1 = a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 >1= a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
: 0.2107209027515327118501721676 * a - 0.2107209027515327118501721676 
: -1.181014228762492651381540270 * a + 1.181014228762492651381540270 
= 2.637456845238095238095238095 * a - 2.637456845238095238095238095 
= 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
-1.667163519227135298563869992 * a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 >1= a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
: 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
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v 8 * W2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 8 *w 3 = 0.09050595238095238095238095238*a— 0.09050595238095238095238095238 
v 8 *w A = — 0.4262670212375073486184597296*a+0.4262670212375073486184597296 
V8 * W5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 8 * We = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 8 * w 7 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
vs *vu 8 = 0.3357610688565549676660787772 * a + 0.0000000000000000000000000000 
V 8 * w 9 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 9 * w x = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 9 * w 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
Vg * w 3 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
vg*w A = 0. 009602757962753002435542118082*0;— 0.009602757962753002435542118082 
Vg * w 5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
Vg * w e = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
Vg * w 7 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
Vg *w 8 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

v 9 *wg = —0. 0096027579627530024355421 18082*a+0. 0000000000000000000000000000 


Flux point j=3 

v\*w\ = — 0.4192916666666666666666666667*a+0. 0000000000000000000000000000 
vi *w 2 = 0.0000000000000000000000000000* a + 0.0000000000000000000000000000 
V! *w 3 = 0.0000000000000000000000000000* a + 0.0000000000000000000000000000 
V\ * w A = 0.8402500000000000000000000000 * a - 0.8402500000000000000000000000 
v\*w§ = — 0.5449166666666666666666666667*a+0. 5449166666666666666666666667 
V! * w 6 = 0.1239583333333333333333333333 * a - 0.1239583333333333333333333333 
V! *w 7 = 0.0000000000000000000000000000* a + 0.0000000000000000000000000000 
V! *w 8 = 0.0000000000000000000000000000* a + 0.0000000000000000000000000000 
v 2 * w x = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 2 *w 2 = 1.458512047141387419165196943 * a + 0.0000000000000000000000000000 
v 2 *w 3 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 2 * w A = -3.422754260886453346770807088 * a + 3.422754260886453346770807088 
v 2 * w 5 = 3.013122649754346182917611489 * a - 3.013122649754346182917611489 
v 2 * w 6 = -1.259601338760812967162173511 * a + 1.259601338760812967162173511 
V2 * W7 = 0.2107209027515327118501721676 * a - 0.2107209027515327118501721676 
V2 * W8 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 3 * Wx = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 3 * w 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 3 *w 3 = —0.5392203804747207524985302763*a+0. 0000000000000000000000000000 
v 3 * w A = 3.746770181143445032333921223 * a - 3.746770181143445032333921223 
v 3 * w 5 = -6.291230553718400940623162845 * a + 6.291230553718400940623162845 
v 3 * w 6 = 4.174189029431216931216931217 * a - 4.174189029431216931216931217 
v 3 * w 7 = -1.181014228762492651381540270 * a + 1.181014228762492651381540270 
v 3 *uj 8 = 0.09050595238095238095238095238*a— 0.09050595238095238095238095238 
v A * Wx = 0.8402500000000000000000000000 * a - 0.8402500000000000000000000000 
v A * w 2 = -3.422754260886453346770807088 * a + 3.422754260886453346770807088 
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v 4 * w 3 = 3.746770181143445032333921223 * a - 3.746770181143445032333921223 
V 4 , * w 4 = - 1 .164265920256991685563114135 * a + 0.0000000000000000000000000000 
v A * w 5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

V A * m = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 4 *w 7 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 4 * w 8 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 3 *w\ = —0.5449166666666666666666666667*a+0. 5449166666666666666666666667 
v 5 * w 2 = 3.013122649754346182917611489 * a - 3.013122649754346182917611489 
v 5 *w 3 = -6.291230553718400940623162845 * a + 6.291230553718400940623162845 

v 5 * u >4 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

v 5 *w 5 = 3.823024570630721424372218023 * a + 0.0000000000000000000000000000 
v 5 * W6 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 5 * w 7 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 
v 5 * w 8 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 6 * Wx = 0.1239583333333333333333333333 * a - 0.1239583333333333333333333333 
v 6 *w 2 = -1.259601338760812967162173511 * a + 1.259601338760812967162173511 
v 6 *w 3 = 4.174189029431216931216931217 * a - 4.174189029431216931216931217 
v 6 * WA = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

V 6 * w 5 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 

V 6 * iu 6 = -3.038546024003737297388091039 * a + 0.0000000000000000000000000000 
v 6 * u, 7 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 6 *w 8 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 7 * Wl = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 7 *w 2 = 0.2107209027515327118501721676 * a - 0.2107209027515327118501721676 
v 7 *w 3 = -1.181014228762492651381540270 * a + 1.181014228762492651381540270 
V7 * WA = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

V 7 * u, 5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

V 7 * u, 6 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

V 7 * W7 = 0.9702933260109599395313681028 * a + 0.0000000000000000000000000000 

v 7 * Wg = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 8 * Wl = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 8 *w 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
v 8 *w 3 = 0.09050595238095238095238095238*a— 0.09050595238095238095238095238 
v 8 *W 4 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 8 * w 5 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V 8 * w 6 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
V8 * w 7 = 0.0000000000000000000000000000 *a + 0.0000000000000000000000000000 

v 8 *w 8 = —0.09050595238095238095238095238*a+0. 0000000000000000000000000000 


Flux point j=2 

vi*wi = —0.07820833333333333333333333333*a+0. 0000000000000000000000000000 

Vi *w 2 = 0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 

vi*w 3 = —0. 3410833333333333333333333333 =t=a+0. 3410833333333333333333333333 
vi * w 4 = 0.8402500000000000000000000000 * a - 0.8402500000000000000000000000 
vi*W5 = — 0.5449166666666666666666666667*a+0. 5449166666666666666666666667 
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V\ * Wq 
V\ * W7 
V 2 * W\ 

V2 * W2 
V 2 *W 3 
V2 * W4 
V2 * W 5 
V 2 * Wq 
V2 * W7 
V 3 *W\ = 
V 3 * W2 
V 3 *W 3 = 
V 3 * W4 
V 3 * Wq 
V 3 * WQ 
V 3 * W7 
V4 * W\ 

V4 * U>2 
V4 * W 3 
V4 * W4 
V4 * Wq 
V4 * Wq 
V4 * W7 
VQ*W\ = 
Vq * W 2 
VQ * W 3 
Vq * W4 
Vq*Wq 
Vq * Wq 
Vq * W7 
Vq * W 1 
Vq * W2 
VQ*W 3 
Vq * W4 
Vq * Wq 
Vq * Wq 
Vq * W7 
V7 * W 1 
V7 * 'Uj'2 
V7 * W 3 
V7 * W4 
V7 * Wq 
V7 * Wq 
V7*W7 = 


0.1239583333333333333333333333 > 1 = a - 0.1239583333333333333333333333 
0.0000000000000000000000000000 > 1 = a + 0.0000000000000000000000000000 
0.0000000000000000000000000000 > 1 = a + 0.0000000000000000000000000000 
0.5782083333333333333333333333 * a + 0.0000000000000000000000000000 
0.8803037138080540858318636096 * a - 0.8803037138080540858318636096 
: -3.422754260886453346770807088 * a + 3.422754260886453346770807088 
: 3.013122649754346182917611489 * a - 3.013122649754346182917611489 
-1.259601338760812967162173511 * a + 1.259601338760812967162173511 
0.2107209027515327118501721676 * a - 0.2107209027515327118501721676 
—0.3410833333333333333333333333*a+0. 3410833333333333333333333333 
0.8803037138080540858318636096 * a - 0.8803037138080540858318636096 
—0.5392203804747207524985302763*a+0. 0000000000000000000000000000 
0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
0.0000000000000000000000000000 > 1 = a + 0.0000000000000000000000000000 
0.0000000000000000000000000000 > 1 = a + 0.0000000000000000000000000000 
0.8402500000000000000000000000 > 1 = a - 0.8402500000000000000000000000 
-3.422754260886453346770807088 * a + 3.422754260886453346770807088 
0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
: 2.582504260886453346770807088 * a + 0.0000000000000000000000000000 
0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
— 0.5449166666666666666666666667*a+0. 5449166666666666666666666667 
: 3.013122649754346182917611489 * a - 3.013122649754346182917611489 
0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
-2.468205983087679516250944822 * a + 0.0000000000000000000000000000 
0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
0.0000000000000000000000000000 > 1 = a + 0.0000000000000000000000000000 
0.1239583333333333333333333333 * a - 0.1239583333333333333333333333 
-1.259601338760812967162173511 * a + 1.259601338760812967162173511 
0.0000000000000000000000000000 > 1 = a + 0.0000000000000000000000000000 
0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
0.0000000000000000000000000000 > 1 = a + 0.0000000000000000000000000000 
: 1.135643005427479633828840178 * a + 0.0000000000000000000000000000 
0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
0.0000000000000000000000000000 > 1 = a + 0.0000000000000000000000000000 
0.2107209027515327118501721676 * a - 0.2107209027515327118501721676 
0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
0.0000000000000000000000000000 * a + 0.0000000000000000000000000000 
0.0000000000000000000000000000 > 1 = a + 0.0000000000000000000000000000 
—0. 21072090275153271 18501721676*a+0. 0000000000000000000000000000 
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Flux point j=l 

v\*w\ = 0.5000000000000000000000000000* a + 0 
Vi *102 = -0.5782083333333333333333333333*a+0 
V!*W3 = — 0.3410833333333333333333333333*a+0 
v 1 * w A = 0.8402500000000000000000000000 * a - 0 
v\*w§ = — 0.5449166666666666666666666667*a+0 
Vi * w 6 = 0.1239583333333333333333333333 * a - 0 
v 2 *wi = — 0.5782083333333333333333333333*a+0 
v 2 *w 2 = 0.5782083333333333333333333333 *a + 0 
v 2 * '(03 = 0.0000000000000000000000000000* a + 0 
v 2 *W4 = 0.0000000000000000000000000000* a + 0 
v 2 *we = 0.0000000000000000000000000000* a + 0 
V2*we = 0.0000000000000000000000000000* a + 0 
V3* Wl = — 0.3410833333333333333333333333*a+0 
v 3 *w 2 = 0.0000000000000000000000000000* a + 0 
v 3 *w 3 = 0.3410833333333333333333333333 *a + 0 
v 3 *W4 = 0.0000000000000000000000000000* a + 0 
v 3 * w 5 = 0.0000000000000000000000000000* a + 0 
v 3 * vje = 0.0000000000000000000000000000* a + 0 
v A * Wl = 0.8402500000000000000000000000 * a - 0 
V4*w 2 = 0.0000000000000000000000000000* a + 0 
V4 * w 3 = 0.0000000000000000000000000000* a + 0 
v 4 *w 4 = — 0.8402500000000000000000000000*a+0 

V4*W5 = 0.0000000000000000000000000000* a + 0 
V 4 * 'We = 0.0000000000000000000000000000* a + 0 

ve*wi = — 0.5449166666666666666666666667*a+0 
ve*w 2 = 0.0000000000000000000000000000* a + 0 
v$*w 3 = 0.0000000000000000000000000000* a + 0 
vg * 11J4 = 0.0000000000000000000000000000* a + 0 
V5 * = 0.5449166666666666666666666667 * a + 0 

v$*we = 0.0000000000000000000000000000* a + 0 
v 6 * Wl = 0.1239583333333333333333333333 * a - 0 
ve*w 2 = 0.0000000000000000000000000000* a + 0 
ve*w 3 = 0.0000000000000000000000000000* a + 0 
ve*W4 = 0.0000000000000000000000000000* a + 0 
ve * 105 = 0.0000000000000000000000000000* a + 0 
v 6 *w 6 = — 0.1239583333333333333333333333*a+0 


.0000000000000000000000000000 

.5782083333333333333333333333 

.3410833333333333333333333333 

.8402500000000000000000000000 

.5449166666666666666666666667 

.1239583333333333333333333333 

.5782083333333333333333333333 

.0000000000000000000000000000 

.0000000000000000000000000000 

.0000000000000000000000000000 

.0000000000000000000000000000 

.0000000000000000000000000000 

.3410833333333333333333333333 

.0000000000000000000000000000 

.0000000000000000000000000000 

.0000000000000000000000000000 

.0000000000000000000000000000 

.0000000000000000000000000000 

.8402500000000000000000000000 

.0000000000000000000000000000 

.0000000000000000000000000000 

.0000000000000000000000000000 

.0000000000000000000000000000 

.0000000000000000000000000000 

.5449166666666666666666666667 

.0000000000000000000000000000 

.0000000000000000000000000000 

.0000000000000000000000000000 

.0000000000000000000000000000 

.0000000000000000000000000000 

.1239583333333333333333333333 

.0000000000000000000000000000 

.0000000000000000000000000000 

.0000000000000000000000000000 

.0000000000000000000000000000 

.0000000000000000000000000000 


Flux point j=0 

Vi * Wl = 0.0000000000000000000000000000 * a + 1.000000000000000000000000000 


The P-norm V is: 

Vij/dx = 0.29950000000000000000000000000000000000 
V 2 ,2 /dx = 1.4888460333207357016880826404635928445 
V 3 , 3 /dx = 0.38651394400352733686067019400352733686 
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V 4A /dx = 1.5399915674603174603174603174603174603 
V^s/dx = 0.73536072530864197530864197530864197531 
Vf- A /dx = 1.0203624889770723104056437389770723104 
V 7 j/dx = 1.0523566468253968253968253968253968254 
Vs.s/dx = 0.97245927028218694885361552028218694885 
V 93 /dx = 1.0046093238221214411690602166792642983 
Vi 0 ,io/dx = 1.0000000000000000000000000000000000000 

The P-norm inverse V~ l is: 

■P u ' * dx = 3.3388981636060100166944908180300500835 

* dx = 0.67166112386355419609800361200830541810 

* dx = 2.5872287805246011319981479868054754454 

* dx = 0.64935420500331280158653723817674480927 

* dx = 1.3598768136281481908969048637106484335 
p-i * dx = 0.98004386754996649194151043751580269422 

* dx = 0.95024819106399044849734934488965855641 

* dx = 1.0283207025316559720897982867102527893 
V 9 g * dx = 0.99541182456421479599641740148968966041 
V w ] w *dx = 1.0000000000000000000000000000000000000 


The Qij * dx nonzero terms are: 

Qi .1 = -0.50000000000000000000000000000000000000 
Qr’ 2 = 0.57820833333333333333333333333333333333 
Qi’ 3 = 0.34108333333333333333333333333333333333 
Q 1)4 = -0.84025000000000000000000000000000000000 
Qi,5 = 0.54491666666666666666666666666666666667 
Qi’ 6 = -0.12395833333333333333333333333333333333 
Q 2 ’r = -0.57820833333333333333333333333333333333 
Q 2 ’ 3 = -0.88030371380805408583186360964138741917 
Q 2)4 = 3.4227542608864533467708070882674057277 
Q 2 ’ 5 = -3.0131226497543461829176114890400604686 
Q 2 ’ 6 = 1.2596013387608129671621735113798605862 
Q 27 = -0.21072090275153271185017216763248509280 
g 3 ’r = -0.34108333333333333333333333333333333333 
Q 3 ’ 2 = 0.88030371380805408583186360964138741917 
Q 3A = -3.7467701811434450323339212228101116990 
Q 3 ’ 5 = 6.2912305537184009406231628453850676073 
g 3 ’ 6 = -4.1741890294312169312169312169312169312 
Q 3J = 1.1810142287624926513815402704291593181 
Q 3 ’ 8 = -0.090505952380952380952380952380952380952 
Gt’i = 0.84025000000000000000000000000000000000 
Q 4 ’ 2 = -3.4227542608864533467708070882674057277 
Q 4 ’ 3 = 3.7467701811434450323339212228101116990 
Q 4 ’ 5 = -5.1042661439961787184009406231628453851 
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Q 4 , 6 = 6.1607928057025279247501469723691945914 
Q 4J = -2.6374568452380952380952380952380952381 
Q 4 ’ 8 = 0.42626702123750734861845972957084068195 
Q 4 ’g = -0.0096027579627530024355421180818006214832 
0 3 ,i = —0.54491666666666666666666666666666666667 
Q 5 ’ 2 = 3.0131226497543461829176114890400604686 
Q 5 ’ 3 = -6.2912305537184009406231628453850676073 
g 5 ’ 4 = 5.1042661439961787184009406231628453851 
Q 5 ’ 6 = -3.3485729166666666666666666666666666667 
Q 5 ’ 7 = 2.7798440990593768371546149323927101705 
Q 5 ’ s = -0.75580278604497354497354497354497354498 
Q 5 ’ 9 = 0.043290030286806080456874107667758461409 
Qe’i = 0.12395833333333333333333333333333333333 
Q 6 ’ 2 = -1.2596013387608129671621735113798605862 
Q 6 ’ 3 = 4.1741890294312169312169312169312169312 
Q 64 = -6.1607928057025279247501469723691945914 
Q 6 ’ 5 = 3.3485729166666666666666666666666666667 
Qtjj = -0.66288489767048794826572604350382128160 
Q 6 ’g = 0.49191676771017048794826572604350382128 
Q 6 ’ 9 = -0.051786576436130007558578987150415721844 
Qe’io = -0.0035714285714285714285714285714285714286 
Q t ’ 2 = 0.21072090275153271185017216763248509280 
Q 7 ’ 3 = -1.1810142287624926513815402704291593181 
Q 7 ,4 = 2.6374568452380952380952380952380952381 
Q 7 ’ 5 = -2.7798440990593768371546149323927101705 
Q 7 ’ 6 = 0.66288489767048794826572604350382128160 
Q 7 ’ 8 = 0.53307226447677836566725455614344503233 
Q 7 ’ 9 = -0.11780039183883429915175946921978668010 
Q 7 ’ 10 = 0.038095238095238095238095238095238095238 
Q 7 ’u = -0.0035714285714285714285714285714285714286 
Q 8 ’ 3 = 0.090505952380952380952380952380952380952 
Q 8 ’ 4 = -0.42626702123750734861845972957084068195 
Qs,5 = 0.75580278604497354497354497354497354498 
Q 8 ’ 6 = -0.49191676771017048794826572604350382128 
0,8,7 = -0.53307226447677836566725455614344503233 
Q 8 ’ 9 = 0.77042350547472075249853027630805408583 
Q 8 ’ 10 = -0.20000000000000000000000000000000000000 
Qs’n = 0.038095238095238095238095238095238095238 
Q 8 ’ 12 = -0.0035714285714285714285714285714285714286 
Q 9 ’ 4 = 0.0096027579627530024355421180818006214832 
Qq,5 = -0.043290030286806080456874107667758461409 
Q 9 ’ 6 = 0.051786576436130007558578987150415721844 
Q 9 ’ 7 = 0.11780039183883429915175946921978668010 
Q 9 ’ 8 = -0.77042350547472075249853027630805408583 
Q 9 ’ 10 = 0.80000000000000000000000000000000000000 


63 



g 9 , n = -0.20000000000000000000000000000000000000 
Q 9 ’ 12 = 0.038095238095238095238095238095238095238 
Q 9 ’ 13 = -0.0035714285714285714285714285714285714286 
QiO, 6 = 0.0035714285714285714285714285714285714286 
Qio’t = -0.038095238095238095238095238095238095238 
Qio’s = 0.20000000000000000000000000000000000000 
Qio’g = -0.80000000000000000000000000000000000000 
Qio’n = 0.80000000000000000000000000000000000000 
Qio .12 = -0.20000000000000000000000000000000000000 
Q 10 ’ 13 = 0.038095238095238095238095238095238095238 
Qio’u = -0.0035714285714285714285714285714285714286 

The (4-8-4) viscous terms. 

The TZi j * dx nonzero terms are: 

K 0 ,i = -2.4500000000000000000000000000000000000 
1z 02 = 6.0000000000000000000000000000000000000 
TZ 0}3 = -7.5000000000000000000000000000000000000 
n 0A = 6.6666666666666666666666666666666666667 
TZo,5 = -3.7500000000000000000000000000000000000 
K 0fi = 1.2000000000000000000000000000000000000 
IZoj = —0.16666666666666666666666666666666666667 
Ki A = -1.1602083333333333333333333333333333333 
7t i i2 = 1.1564166666666666666666666666666666667 
^i’ 3 = 0.34108333333333333333333333333333333333 
= —0.56016666666666666666666666666666666667 
ni i5 = 0.27245833333333333333333333333333333333 
nifi = -0.049583333333333333333333333333333333333 
n 2 , i = 0.093973038486880787332060994325567220421 
K-2 , 2 = -0.78549893930185917599538939754826348881 
TZ 2:3 = 0.046946483020436971650521028934066801314 
1l 2A = 0.90729349114883760012958720175606131931 
TZ2,5 = -0.26981952253301897863082607780986533944 
n 26 = -0.0063708948742115737446125435968063627529 
K 2 ,7 = 0.013476344052934369258658793939239849966 
n 3A = 0.078746435191005348453800791189557639584 
77.3/2 = -0.35162113760879514268346232087309277418 
1l 33 = -0.78713436791133549855006393774752877572 
U 3A = 1.3893360044549165552850116426510172239 
77 3 ’ 5 = -0.45575020893114105530342782230270280369 
1l 3fi = 0.19808018162235668097172402106009418736 
TZ 3J = -0.092080135398545507498723070166851812875 
77 3 ’ 8 = 0.020423228581538619325140696189507115568 
TZ 4 a = 0.032926148937526826310667271834304379849 
K A2 = -0.25191802597411192851286613885119855940 
= 0.86404038300564659933759665503107585561 
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TZ 4 a = -2.2705734372043358134175184621100160529 
n 4 ,5 = 1.9056036402801850141544340055650364152 
H 46 = -0.43544242414791128227318211540346767698 
n 4 , 7 = 0.20103197771146276694189556489629198797 
ft 4 ’ 8 = -0.049509365793563383515243628194746597938 
= 0.0038411031851012009742168472327202485933 
K bA = -0.049583333333333333333333333333333333333 
77.5/2 = 0.24538547964509061280525769040437984250 
= -0.46041931266967499182355759462391059263 
77 . 5 ’zt = 0.46614273043637761869004215693100311426 
77 5 ’ 5 = -1.3962385850218969553701412457089944128 
77 5 /j = 1.3149492005817200178282234837505867257 
77,5 j = -0.13710317505277669221310366487172969582 
77 5 ’ 8 = 0.025083621435175386943099478252132619196 
ft 5 ’g = -0.0082166260206816635264869708001342671278 
77 6 ’ 2 = 0.013476344052934369258658793939239849966 
77f/ 3 = -0.095466650296571124068030350198485565808 
77 6 ’ 4 = 0.28303601526589854838736615694592248492 
77 6 ’ 5 = -0.37189053652589385913924939420370302618 
77 64 j = -0.87133674643234682661783233764660965659 
77 6 ’ 7 = 1.1363503084863723276987369694549351287 
^ 6 ’ 8 = -0.10962852877844960704052530455432144803 
77. 6 /9 = 0.017245508513770457235161180548736518789 
77 6 ’i 0 = -0.0017857142857142857142857142857142857143 
77 . 7^3 = 0.020423228581538619325140696189507115568 
77.7^4 = -0.11461749320859847500014290897545707151 
ft.7’5 = 0.24553622332394432725008432375299558042 
77.7/j = -0.14002362171685620274536032506670231808 
ft.7’7 = -1.2086425337884685391310360256213220479 
n 7 [ 8 = 1.3257988165251715762302257817106939321 
n 7 ,g = -0.15030001654212813132573693881511201603 
n 7A0 = 0.023611111111111111111111111111111111111 
U 7M = -0.0017857142857142857142857142857142857143 
K 8 ,4 = 0.0038411031851012009742168472327202485933 
1Z 8t5 = -0.021645015143403040228437053833879230705 
1l 8fi = 0.036310098576467624086671705719324766944 
n 8J = 0.094189280727723188040648358108675568993 
1Z 8 .s = -1.3644581220605526161081716637272192828 
1l 8:9 = 1.4063261467781557067271352985638699924 
TZ 8 ao = -0.17638888888888888888888888888888888889 
Ks.n = 0.023611111111111111111111111111111111111 
U 8A2 = -0.0017857142857142857142857142857142857143 
77-9 /j = 0.0017857142857142857142857142857142857143 
n 97 = -0.023611111111111111111111111111111111111 
Hg :8 = 0.17638888888888888888888888888888888889 
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7l 9 ,9 = -1.4236111111111111111111111111111111111 
ftg’io = 1.4236111111111111111111111111111111111 
ftg’n = -0.17638888888888888888888888888888888889 
H 912 = 0.023611111111111111111111111111111111111 
H 9t i 3 = -0.0017857142857142857142857142857142857143 


The TZ h j * dx nonzero terms are: 

TZ h i = -1.1602083333333333333333333333333333333 
TZi t 2 = 1.1564166666666666666666666666666666667 
lZi£ = 0.34108333333333333333333333333333333333 
1Z 4 , 4 = —0.56016666666666666666666666666666666667 
Tlifi = 0.27245833333333333333333333333333333333 
nifi = -0.049583333333333333333333333333333333333 
n-2,i = 1.2541813718202141206653943276589005538 
712,2 = -1.9419156059685258426620560642149301555 
K 2 , 3 = -0.29413685031289636168281230439926653202 
K 2 a = 1.4674601578155042667962538684227279860 
n 2 ,5 = -0.54227785586635231196415941114319867278 
= 0.043212438459121759588720789736526970580 
n 2 , 7 = 0.013476344052934369258658793939239849966 
TZ 3 ,i = -0.015226603295875438878260203136009580837 
n 3:2 = 0.43387780169306403331192707667517071464 
U 3 , 3 = -0.83408085093177247020058496668159557703 
n 3A = 0.48204251330607895515542444089495590463 
7Z 3:5 = -0.18593068639812207667260174449283746424 
U 3fi = 0.20445107649656825471633656465690055011 
n 3J = -0.10555647945147987675738186410609166284 
U 3 ,s = 0.020423228581538619325140696189507115568 
= -0.045820286253478522143133519355253259735 
nit = 0.099703111634683214170596182021894214773 
ni r 3 = 1.6511747509169820978876605927786046313 
U 4 , 4 = -3.6599094416592523687025301047610332768 
K 4 , 5 = 2.3613538492113260694578618278677392189 
= -0.63352260577026796324490613646356186434 
n 4J = 0.29311211311000827444061863506314380084 
n 4 , s = -0.069932594375102002840384324384253713506 
n 4 < 4 = 0.0038411031851012009742168472327202485933 
= -0.082509482270860159644000605167637713182 
n 4 - 2 = 0.49730350561920254131812382925557840190 
TZ b:3 = -1.3244596956753215911611542496549864482 
n 5 , 4 = 2.7367161676407134321075606190410191672 
U 5 ’ 5 = -3.3018422253020819695245752512740308280 
= 1.7503916247296313001014055991540544027 
U b j = -0.33813515276423945915499922976802168379 
Ko.s = 0.074592987228738770458343106446879217134 
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U bt 9 = -0.012057729205782864500703818032854515721 
K 6 , i = 0.049583333333333333333333333333333333333 
7e 6 ’ 2 = -0.23190913559215624354659889646513999253 
U 6:3 = 0.36495266237310386775552724442542502682 
Kf- A = -0.18310671517047907030267599998508062934 
= 1.0243480484960030962308918515052913866 
= -2.1862859470140668444460558213971963823 
= 1.2734534835391490199118406343266648245 
7l 6: 8 = -0.13471215021362499398362478280645406723 
= 0.025462134534452120761648151348870785916 
ft 6 ’ 10 = -0.0017857142857142857142857142857142857143 
n 7 [ 2 = -0.013476344052934369258658793939239849966 
U 7 [ 3 = 0.11588987887810974339317104638799268138 
n 7 . A = -0.39765350847449702338750906592137955643 
U 7 [ 5 = 0.61742675984983818638933371795669860660 
U 7fi = 0.73131312471549062387247201257990733851 
n 7 [ 7 = -2.3449928422748408668297729950762571765 
TZ 7 [ s = 1.4354273453036211832707510862650153801 
H 7 $ = -0.16754552505589858856089811936384853482 
U 7 ] io = 0.025396825396825396825396825396825396825 
n 7U = -0.0017857142857142857142857142857142857143 
TZ 8}3 = -0.020423228581538619325140696189507115568 
U 8A = 0.11845859639369967597435975620817732010 
1Z 8A = -0.26718123846734736747852137758687481113 
U 8fi = 0.17633372029332382683203203078602708503 
U 8 j = 1.3028318145161917271716843837299976169 
TZ 8t8 = -2.6902569385857241923383974454379132149 
K 8A = 1.5566261633202838380528722373789820085 

n 8A o = -0.20000000000000000000000000000000000000 

U 8M = 0.025396825396825396825396825396825396825 
n 8A2 = -0.0017857142857142857142857142857142857143 
K 9A = -0.0038411031851012009742168472327202485933 
Ug A = 0.021645015143403040228437053833879230705 
U 9fi = -0.034524384290753338372385991433610481230 
Kgj = -0.11780039183883429915175946921978668010 
Kg, 8 = 1.5408470109494415049970605526161081717 
U 9A = -2.8299372578892668178382464096749811035 
K 9A0 = 1.6000000000000000000000000000000000000 
ftg’n = -0.20000000000000000000000000000000000000 
n 9A2 = 0.025396825396825396825396825396825396825 
U 9A3 = -0.0017857142857142857142857142857142857143 
TZ 10fi = -0.0017857142857142857142857142857142857143 
K 10 j = 0.025396825396825396825396825396825396825 
= -0.20000000000000000000000000000000000000 
ftio’g = 1.6000000000000000000000000000000000000 
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Ki 0 ,io = - 2.8472222222222222222222222222222222222 
TZ w ,n = 1.6000000000000000000000000000000000000 

■Rio t i2 = - 0.20000000000000000000000000000000000000 

n 10 ,!3 = 0.025396825396825396825396825396825396825 
n 10M = - 0.0017857142857142857142857142857142857143 


The BSij nonzero terms are: 

BS hl = 2.4500000000000000000000000000000000000 
BS h2 = - 6.0000000000000000000000000000000000000 
BS h 3 = 7.5000000000000000000000000000000000000 
BSi^ = — 6.6666666666666666666666666666666666667 
BSi )5 = 3.7500000000000000000000000000000000000 
BSifi = - 1.2000000000000000000000000000000000000 
BSij = 0.16666666666666666666666666666666666667 
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